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Abstract 

Large-scale velocity fields in the solar photosphere remain a mystery in spite of many years of 
intensive studies. In this thesis, the new method of the measurements of the solar photospheric 
flow fields is proposed. It is based on local correlation tracking algorithm applied to full-disc 
dopplergrams obtained by Michelson Doppler Images (MDI) on-board the Solar and Heliospheric 
Observatory (SoHO). The method is tuned and tested on synthetic data, it is shown that the 
method is capable of measuring of horizontal velocity fields with an accuracy of 15 ms^^. It 
is also shown that the method provides the measurements comparable with time-distance local 
helioseismology. The method is applied to real data sets. It reproduces well known properties 
of solar photospheric velocity fields. Moreover, the case studies show an evidence about the 
influence of the changes in the flow fleld topology on the stability of the eruptive fllament and 
support the theory of the dynamical disconnection of bipolar sunspots from their magnetic roots. 
The method has a great perspective in the future use. The meridional flux transportation speed 
is also studied and it is shown that the direct measurement may differ from time-distance local 
helioseimology in the areas occupied by the strong magnetic field. This result has an impact to 
the flux transport dynamo models, which use the meridional speed as the essential observational 
input parameter. 



Mottoe.'^: 



If at first you don't succeed, try, try, try, try, try, try, try, try again. 

Jack O'Neill, Stargate SG-1, Window of Opportunity 

If you don't follow your dreams, you might as well he a vegetable. 

Burt Munro, The World's Fastest Indian 
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1 Preface 

This thesis presents the collection of results obtained during my Ph.D. study program on 
the topic of surface velocity fields in the photosphere of the Sun and connected topics. 
Rather than the topic as a whole, some subtopics were studied. The main result of this 
work is a data set, which can and will be used for other analyses during next years. 

In Section [2] I present a short overview about the physics and known properties of the 



solar surface flow fields. In Section 2.2.1 the principle of the local correlation tracking 
method, which was used in the most of the work done, is explained. The method was 
applied to full-disc dopplergrams. It was tuned using synthetic data. For this purpose, 
a code generating synthetic solar dopplergrams with known properties was developed. It 
is described in Section [3] as well as the optimal values of free parameters involved in the 
method, and the noise and accuracy estimation. 

The performance of the method is verified comparing its results with results of the 
modern-most method of measurements of surface and sub-surface flow fields. In Section H] 
it is shown that both methods applied to the same set of data provide comparable results. 
The processed dataset and the visualization of the results used during the analysis is 
described in Section [H 

In the next Sections, the particular results obtained with the proposed method are 
described: The long-term properties and periodicities in Section |6] and the flows under the 
eruptive filament and their evolution during the eruptive phase in Section [7j 

In the Section |8| a different but very similar topic is studied. The meridional flux 
transport process is an essential property of the solar global dynamo performance. Using 
a different dataset (magnetic butterfly diagram) and a different technique, the meridional 
flux transportation speed is measured and compared with the results of the time-distance 
helioseismology. It is shown that the surface measurements are biased by motions around 
local magnetic regions. This result is important for the flux transport dynamo models. 

The flow fields in the previous Sections concerned only the surface of the Sun. In 



Section 10 it is shown that the local correlation tracking method could work also on stellar 
data coming from the Doppler imaging. So far, the study of this problem is in the very 
early stage, therefore it is attached as an appendix. However, it has a large perspective 
into the future because of lots of data coming from robotic telescopes. 
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2 Introduction 

The Sun is the closest star - this fact allows us to resolve individual features on its surface 
and in its atmosphere. Using many types of observations, we can collect a large amount 
of data describing the behaviour of the solar plasma in various phenomena. The Sun is a 
variable star - the magnetic activity undergoes the main cycle with a length of 22 years 
(reversal of a global magnetic field) which is composed of two consecutive 11 years 'spot' 
cycles. The most visible evidence of the solar cyclicity is the change of the number, size, 
and shape of sunspots. However, evidences of such cyclicity may be found also in the total 
solar irradiance, number of solar flares, or the shape of the solar corona. The Sun exhibits 
also cycles with different lengths (from few minutes to many centuries) and different prop- 
erties. Better knowledge of the physics lying under solar magnetic variability and active 
phenomena will improve our attempts to predict solar activity. 



2.1 Structure and dynamics of the photospheric velocity field, espe- 
cially in large-scale 

The first evidence about the velocity fields in the solar photosphere comes from Christoph 
Scheiner, who in 1630 noticed that sunspots near the equator traverse the solar disc faster 
than sunspots in higher latitudes. Carrington (1859) used series of sunspot drawings 
to infer the differential rotation rate and the inclination of the solar rotation axis. Since 
Carrington's measurement of the differential rotation this phenomenon was confirmed many 



times using many techniques. For details see an older review by Schroter ( 1985 ) or a recent 



review by Beck (2000). The existence of the differential rotation was the first evidence 



about the movements of objects in the solar photosphere. Since then many other types of 
motions were detected. 

The energy coming from a thermonuclear fusion in the solar core is from approx. 0.7 i?© 
carried by the convection. The dynamic behaviour in the solar photosphere is therefore 
mainly driven by the plasma motions in the underlying convection zone. The most evi- 
dent manifestation of the convective behaviour is the solar granulation. It is observed in 
white-light as a cellular pattern with a characteristic size of 1,000 km and a lifetime of 
3-17 minutes (Dialetis et al. , 1986). Motions of granules are studied mainly using local 



correlation tracking technique; see e. g. recent study using the TRACE white light data by 



Krijger et al. (2002). From such studies it is known that granules are carried by the flow 



field of the larger scale and form the cellular-like pattern of supergranulation. On smaller 
scales, exploding granules seem to form cells smaller than supergranules - the mesogran- 



ulation - that are also advected within the supergranular cells (Leitzinger et al. , 2005). 



Nevertheless, the convection power spectrum in the photosphere obtained from Doppler 
measurements peaks at granulation scales il > 1000, I is a spherical harmonic degree), 
with a secondary peak at / ~ 120 corresponding to supergranulation. There is no evidence 



(Hathaway et al. , 2000) about the peak corresponding to mesogranules or other cellular- 



like convection mode. From this reason, the existence of mesogranulation as a separated 
convection mode is still in doubt. 
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2.1.1 Supergranules 



The supergranulation pattern is clearly seen as a peak in a convection power spectrum with 



/ ~ 120 (Hathaway et al. 2000). This corresponds to size of roughly 30 Mm. Supergranules 



were discovered by Hart (1956) when studying variations of the equatorial rotation velocity 



using the autocorrelation method. He found a pattern with a characteristic size of 26 Mm 
and a characteristic velocity amplitude of 170 m s~^. A more detailed study done using full- 
disc dopplergrams was performed by Leighton et al. (1962), who basically confirmed Hart s 



results. In the continuing studies by Noyes & Leighton ( 1963 ) and Simon & Leighton ( 1964 ) 



also for the first time appeared the name of the new convection mode - "super-granulation" . 
After a decades of the studies, the true origin of supergranulation remains a mystery. There 
exist many realistic simulations of the solar convection zone, however in none of them it was 
reported the convection mode with the properties of supergranules. Convective motions 



at supergranular scales have been reported in global simulations by DeRosa et al. (2002), 



who focused on the upper regions of the solar convection zone. Higher spatial resolution 
than in other global simulations was achieved by limiting the simulation domain to radii 
between 0.92-0.98 Rq and by imposing a four-fold periodicity in longitude. The resulting 
pattern exhibits a hierarchy of scales, from supergranular-scale mottling to a network of 
larger cells and extended north-south downflow lanes more comparable to the deep-shell 
simulations. Although provocative, it is premature to identify this small-scale convection 
pattern too closely with supergranulation on the Sun, because other statistical parametres 
of the results of simulation do not agree with the parametres of real supergranulation. Solar 
supergranulation may involve dynamics, which is not captured in these global simulations 
such as ionization effects or self-organization processes involving smaller-scale granules 



(Rast, 2003). On the other hand, there is no evidence of supergranulation in the detailed 



simulations focused on granulation. The origin of supergranulation remains unclear, Rincon 



et al. (2005) proposed the collective interaction of the granules, as its cause. In particular. 



the author showed that the formation of the long-lived large scale pattern can be obtained 
by computing the advection of many small-scale short-lived granular downfiows. 

The determined size of the supergranular cells is sensitive to a method used for the 
measurement of this quantity. Most often it is used the value of 30 Mm. Hart (1956), 
using the autocorrelation method, got the size of 26 Mm. Wang & Zirin (1989) got in 



a detailed study using the autocorrelation method of full-disc dopplergrams the size of 



31.2±2.3 Mm. Srikanth et al. (2000) used a tesselation method to get a typical size of 
the supergranule. They applied the method to full-disc MDI dopplergrams (2"px^^) and 
also to Ca II K filtergrams obtained at the South Pole Observatory (3.2" px^^). They 
found 10.5 Mm in the photosphere and 14-26 Mm in the chromosphere as a typical size 
of the supergranule. However, they argued that the results obtained by the tesselation 
method are strongly influenced by the resolution of images. Using statistical tests they 
concluded that the most probable size of the supergranule in the photosphere and the 
chromosphere is 25.9 Mm. Also Hagenaar et al. (1996) noted that the autocorrelation 



method has a tendency to overestimate (1.5-2 times) the size of real structures. The depth 
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of the supergranulation, inferred using local helioseismology and using various studies, is 



supposed to have a value of 8 Mm (Duvall, 1998) to 15 Mm (Zhao & Kosovichev, 2003). 



The lifetime of supergranulation is one of the crucial parameters. One of the early 
studies (Wang & Zirin, 1989) used the dopplergrams measured at NSO Kitt Peak. They 



segmented the series of images, labeled the supergranules and tracked each of them in 
the whole series. Using exponential fit they obtained the characteristic lifetime of 50-80 
hours. They also mentioned some of cells that survived over one week. In the recent study 



by DeRosa et al. (2000), which is a part of DeRosa's thesis, the authors made a carefull 



study of supergranular lifetimes, made a histogram of their lifetimes, and reported that 
the long-duration datasets contain several instances, where individual supergranules are 
recognizable for time scales as long as 50 hours, though most cells persist for about 25 



hours that is often quoted as a supergranular lifetime. Del Moro et al. (2004) found the 



mean lifetime of 22.5 hours. They also studied the relation between the lifetime and the 
size of the supergranules and found that supergranules probably exist in two regimes. For 
small supergranules (under 27 Mm) there can be observed almost linear increase of the 
mean lifetime with size, for larger (above 27 Mm) is the mean lifetime nearly constant 
(33 hours). This behaviour can be related to the loss of structural coherence by the largest 
supergranules, leading to a fragmentation event in smaller parts. 

The internal velocity field in supergranules is nearly horizontal and is not easy to be 



measured. E.g. Hathaway et al. (2002) used full-disc MDI dopplergrams with solar rota- 



tion, meridional circulation and p-modes of solar oscillations removed. They decomposed 
the dopplergram in the series of annuli of different radius r, which is related to the helio- 
centric angle p = arcsin r/i?, where R is the solar radius. Authors assumed that the mean 
line-of-sight velocity Vd in the annulus may be calculated using 



+ 



- in 



sin^ p, 



(1) 



where v^j is the vertical and Vh the horizontal component of the internal supergranular 
velocity field. Using the least-square fit they obtained the typical values of the internal 
horizontal velocity of (vh) = (258 ± 1) ms~^ and of the internal vertical velocity of (vy) = 



(29 ± 2) ms . It is coherent with many older studies (e. g. Leighton et al. , 1962) 



Simon & Leighton (1964) noted that supergranule boundaries correspond with the 



chromospherical Ca II K network. This behaviour is coherent with the idea of the advected 
magnetic flux tubes, which are advected towards the boundaries of supergranular cells and 
concentrate there. Similar network can be observed also in other 'metal' lines. 



2.1.2 Giant cells? 

There is no doubt that large-scale velocity structures with I < 100 exist in the convection 
power spectra obtained from surface measurements or local helioseismology. However, there 
has not been found a clear evidence about the characteristic pattern of such structures or 
obtained the qualitative diagnostics of large-scale convective motions. In the velocity power 
spectra, there are two peaks denoting existence of granulation and supergranulation. At 
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low / the power of velocity mode at / drops almost linearly with / (Hathaway et al. , 2000). 
Of course this fact could also mean that there is no characteristic pattern in large-scale. 



The giant cells were noted is some older studies (Bumba, 1970 or Bumba, 1987) con- 



nected with the distribution of the magnetic field in the solar photosphere. Also Ambroz 



(1997) used low-resolution magnetic maps as tracers and found some evidence about the ex- 
istence of giant cells. Moore et al. (2000) found pattern with characteristic size 3-10-times 
larger than the size of supergranules with lifetime greater than 10 days. When analysing 
power spectra of full-disc MDI dopplergrams, they concluded that the physical origin of 
such giant cells and the supergranulation is the same - 'their' giant cells are just larger 
and more stable. They also found the correspondence between 'large-scale supergranules' 
and the length of large and stable bipolar sunspot groups. 

Recently, several groups have reported long-lived features in dopplergrams which are 
highly correlated in longitude, corresponding to azimuthal wavenumbers of m = 0-8 (an- 
gular extent > 45°) but with a narrow latitudinal extent of not more than about 6°. 
Although Beck et al. (1998) interpret these features as giant convection cells, Ulrich (2001) 
argues that they more likely comprise a spectrum of inertial oscillations, possibly related 
to Rossby wave modes and perhaps also to torsional oscillations. Evidence for a different 
giant cell structure has been presented by Lisle et al. (2004). They studied the supergranu- 
lation pattern using correlation tracking and found a tendency for north-south alignment of 
supergranular cells. Such an alignment would be expected if the supergranulation were ad- 
vected by larger-scale, latitudinally-elongated lanes of horizontal convergence such as those 
commonly seen in numerical simulations of solar convection (Brun et al. , 2004). Advection 
by such structures may also help to explain why the supergranulation pattern appears to 
rotate faster than the surrounding plasma measured by Doppler shift. 

We can conclude that there is still not any clear evidence about the existence of the 
giant cells. It is assumed that the size of such cells should be between 200 and 400 Mm, 
lifetime greater than one week and the internal velocity field is expected mostly horizontal 
with the magnitude of few meters per second. Early numerical simulations (Simon & 



Weiss, 1991) showed that the existence of giant cell convection is not neccessary for the 
heat transport within the whole convection zone. They concluded that the existence of 
convection cells with size of 10-20 Mm forming at the base of the convection zone with 
a distance of 30-40 Mm between neighbouring cells is sufficient for heat transport to the 
surface. From such travelling cells, the supergranulation is formed in the shallow subsurface 



layer. Recent global simulations by Brun et al. (e. g. 2004) show global convection of many 



kinds, parameters of whose are sensitive to the state parameters of plasma in simulations. 

Global scale convection is observed using interferometric techniques in the close super- 
giant star Betelgeuse (a Orionis) - see Buscher et al. (1990). 



2.1.3 DifFerential rotation 



Differential rotation is an obvious, yet poorly understood, solar phenomenon. Systematic 
measurements of its parametres were started by Carrington (e.g. 1859), who derived the 
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regression formula describing the differentiality of the spot rotation in form 

uj = A + B sin^/^ 6, (2) 

where b is the hehographic latitude and uj is the angular rotation rate. 

Since then, there exist many methods of measuring the solar rotation. They have one 
common denominator - the results differ with the dataset, with the method, and with the 
resolution. 

Basically, the differential rotation is described as an integral of the zonal component 
of the studied flow field. The integrated flow field may be obtained using spectroscopic 
method, using tracer-type measurements, or using helioseismic inversions. The latest case 
allows to measure the solar rotation not only as a function of heliographic latitude, but 
also as a function of depth. From the helioseismic inversion we know that throughout the 
convective envelope, the rotation rate decreases monotonically toward the poles by about 
30 %. Angular velocity contours at mid-latitudes are nearly radial. Near the base of the 
convection zone, there is a sharp transition between differential rotation in the convective 
envelope and nearly uniform rotation in the radiative interior. This transition region 
has become known as the solar tachocline. The rotation rate of the radiative interior 
is intermediate between the equatorial and polar regions of the convection zone. Thus, 
the radial angular velocity gradient across the tachocline is positive at low latitudes and 



negative at high latitudes, crossing zero at a latitude of about 35 ° (e. g. Thompson et al. 



2003). In addition to the tachocline, there is another layer of comparatively large radial 
shear in the angular velocity near the top of the convection zone. At low and mid-latitudes 
there is an increase in the rotation rate immediately below the photosphere, which persists 
down to r ~ 0.95 Rq. The angular velocity variation across this layer is roughly 3 % of 



the mean rotation rate and according to the helioseismic analysis of Corbard &: Thompson 



(2002) a; decreases within this layer approximately as r~^. At higher latitudes the situation 
is less clear. The radial angular velocity gradient in the subsurface shear layer appears to 
be smaller and may switch sign. 

Currently, the surface measurements of the solar rotation are expressed in the form 

u = A + Bsin^b + Csm^b . (3) 

In this formula, A is equivalent to the equatorial rotation rate, B and C describe the 
differentiality. The difficulty of such expression lies in the coupling of B and C by an 
inverse correlation. It makes some issues when comparing different measurements made 
using different methods and different datasets. One remedy to this problem is to fix the 
ratio B/C. However, this technique is not satisfactory enough for comparison of very 
different types of measurements. Another approach, which completely solves the problem, 
is to use as a base function of the fit not even powers of sin b, but to use an orthogonal 



base functions, such as Gegenbauer polynomials (Snodgrass, 1984). 

Generally, the measured solar rotation is more rigid when measured using larger-scale 
objects, such as coronal holes or large-scale background magnetic field. There are known 
many relations of the differential rotation profile to the phase of the progressing solar cycle 
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- see e.g. Javaraiah (2003). For overview of the solar differential rotation measurements 
see 



Schroter (1985) or a more recent review by Beck (2000). 



Current (magneto)hydrodynamic simulation basically reveal the structure of the solar 
rotation. They use many hydrodynamic quantities to reproduce the measured rotation 
profile - angular momentum transport by convection, Reynolds stresses, and also angular 
momentum transport by meridional circulation, which seems to be essential for many 
types of simulations. See e.g. Kiiker et al. (1993) or Brun et al. (2004). However, the 
results are not so satisfying. On the positive side, the angular velocity exhibits a realistic 
latitudinal variation, with little radial variation above mid-latitudes. On the negative side, 
the low-latitude angular velocity contours are somewhat more cylindrical than suggested 
by helioseismology, with more radial shear. Furthermore, at present there is little tendency 
for simulations such as these to form rotational shear layers near the top and bottom of 
the convection zone. Although these simulations do exhibit non-periodic angular velocity 
fluctuations of about the right amplitude relative to helioseismic inversions (a few percent; 
see Miesch, 2000), there is currently little evidence for systematic behavior such as torsional 



oscillations. Since the radiative interior possesses much more mechanical and thermal 
inertia than the convective envelope, the differential rotation in the convection zone may 
be sensitive to the complex dynamics occuring in the tachocline. In other words, we 
may not fully understand the rotation profile in the convection zone until we get the 
tachocline right. A realistic tachocline is probably also a prerequisite to achieving the 
solar-like dynamo cycles and wave-mean flow interactions which appear to be responsible 
for torsional and tachocline oscillations. 



2.1.4 Oscillations in the rotation pattern 



The "torsional oscillations" , in which narrow bands of faster than average rotation, inter- 
preted as zonal flows, migrate towards the solar equator during the sunspot cycle, were 



discovered by Howard & Labonte (1980). At latitudes below about 40°, the bands propa- 



gate equatorward, but at higher latitudes they propagate poleward. The low-latitude bands 
are about 15 ° wide in latitude. The flows were studied in surface Doppler measurements 



(Ulrich, 2001), and also using local helioseismology (Kosovichev & Schou, 1997). The sur- 



face pattern of torsional oscillations penetrate deep in the convection zone, possibly to 
its base, as studied by Vorontsov et al. (2002). The magnitude of the angular velocity 
variation is about 2-5 nHz, which is roughly 1 % (5-10 ms~^) of the mean rotation rate. 
The direct comparison between different techniques inferring the surface zonal flow pattern 



(Howe et al. , 2006) showed that the results are pretty coherent. 



The surface magnetic activity corresponds well with the oscillation pattern - the mag- 
netic activity belt tend to lie on the poleward side of the faster-rotating low-latitude bands. 
The magnetic activity migrate towards the equator with the low-latitude bands of the tor- 



sional oscillations as the sunspot cycle progresses (Zhao h Kosovichev, 2004). Some studies 



e.g. Beck et al. , 2002) suggest that meridional flows may diverge out from the activity 



belts, with the equatorward and poleward flows well correlated with the faster and slower 
bands of torsional oscillations. 
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The other type of oscillations detected in the rotation pattern were first reported by 



Howe et al. (2000). It is completely different from the torsional oscillations. It has a period 



of 1.3 years and the origin of them is localized around the tachocline at the base of the 
convection zone. This type of oscillations depict as a periodic change of the rotation rate, 
the amplitude is about 3 nHz at the equator and slightly larger at a latitude of 60°. So 
far, there is no evidence of latitudinal propagation. The study of oscillations near the 
tachocline is currently at the limits of the sensitivity of helioseismic inversions. However, 



Toomre et al. (2003) found evidence about the 1.3 yr oscillations with variations of the 



amplitude as the solar cycle progresses. 



2.1.5 Meridional circulation 

The differential rotation is the axisymmetric component of the mean longitudinal flow, 
< >. The axisymmetric flow in the meridional plane, < vg > and < Vr >, is generally 
known as the meridional circulation. The meridional circulation in the solar envelope 
is much weaker than the differential rotation, making it relatively difficult to measure. 
Although it can in principle be probed using global helioseismology, the effect of meridional 
circulation on global acoustic oscillations is small and may be difficult to distinguish from 
rotational or other effects. 

Two principal methods are used to measure the meridional flow: feature tracking and 
direct Doppler measurement. There are several difficulties complicating the measurements 
of the meridional flow using tracers. Sunspots and filaments do not provide sufficient 
temporal and spatial resolution for such studies. Sunspots also cover just low latitudinal 
belts and do not provide any informations about the flow in higher latitudes. Doppler 
measurements do not suffer from the problem associated to the tracer-type measurements, 
however they introduce another type of noisy phenomena into account. It is difficult to 
separate the meridional flow signal from the variation of the Doppler velocity from the disc 
center to the limb. Using different techniques the parametres of the meridional flow show 
large discrepancies. It is generaly assumed that the solar meridional flow is in the close 
subphotospherical layers poleward with one cell. Such flow is also produced by early global 



hydrodynamical simulation such as Glatzmaier & Gilman (1982). As reviewed by Hath- 



away (1996), the surface or near sub-surface velocities of the meridional flow are generally 
in range 1-100 ms~^, the most often measured values lie in the range of 10-20 ms~^. The 
flow has often a complex latitudinal structure with both poleward and equatorward flows. 



multiple cells, and large asymmetries about the equator. Zhao & Kosovichev (2004) used 
the time-distance helioseismology to infer the properties of the meridional flow in years 
1996-2002. They found the meridional flows of an order of 20 ms~^, which remained pole- 
ward during the whole period of observations. In addition to the poleward meridional flows 
observed at the solar minimum, extra meridional circulation cells of flow converging toward 
the activity belts are found in both hemispheres, which may imply plasma downdrafts in 
the activity belts. These converging flow cells migrate toward the solar equator together 



with the activity belts as the solar cycle evolves. Beck et al. (2002) measured the merid 



ional flow (and torsional oscillations) using the time-distance helioseismology and found 
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the residual meridional flow showing divergent flow patterns around the solar activity belts 
below a depth of 18 Mm or so. It needs to be noted that the reverse flow (equatorward) , 
which is assumed to flow at the base of the convection zone with velocities of ~ 1 ms~^ 
was not observed yet, although e. g. Hathaway et al. (2003) interpret the migration of the 
sunspot groups during the sunspot cycle as evidence of the deep recurrent meridional flow 
with the average magnitude of 1.2 ms~^. 

Recent global simulations (e.g. Brun et al. , 2004) also reveal the pattern of the merid- 



ional flow. Using the input parametres variation, the agreement between the model and 
measured velocities can be obtained. The modern dynamo flux-transport models use the 
meridional flow and the differential rotation as the observational input. In the models by 
Dikpati et al. (iDikpati et all |2006cl iDikpati k Gilmanl 120061 or iDikpati et all |2006bl) 



the 'retrograde' meridional flow at the base of the convection zone is calculated from the 
continuity equation. They found the turnover time of the single meridional cell of 17- 
21 years. The meridional flow is assumed to be essential for the dynamo action, global 
magnetic fleld reversal and forecast of the ongoing solar cycles. Using the measurements 
of the sub-surface meridional flow they are able to reconstruct the magnetic activity in the 
past cycles and also predict the activity in the ongoing cycles. 



2.1.6 Solar subsurface weather 



The most substantial recent advance in the search for large-scale non-axisymmetric motions 
in the solar envelope has been the mapping of horizontal flows by local helioseismology. 
After subtracting the contributions from differential rotation and meridional circulation, 
the residual flow maps reveal intricate, evolving flows on a range of spatial scales (e. g. Zhao 



& Kosovichev 




2004 


)■ 


(SSW; 


Toomre 


2002) 



Such flow patterns have become known as solar subsurface weather 
(SSW; Toomre, 2002). Such maps derived using different local helioseismic inversions (ring 



The inferred SSW patterns show a high correlation with magnetic activity, becoming 
more complex at solar maximum. Near the surface, strong horizontal flows converge into 
active regions and swirl around them, generally in a cyclonic sense (counter-clockwise in 
the northern hemisphere and clockwise in the southern hemisphere). Deeper down, roughly 
10 Mm below the photosphere, the pattern reverses; here flows tend to diverge away from 
active regions. The topology of the measured SSW is poorly reproduced in recent global 
simulations. The knowledge of the long-term behaviour of the SSW is essential in the fleld 
of investigation of the coupling between velocity and magnetic fleld and contribute in the 
theory of the solar dynamo. 

The inferrence of the 'surface weather' using surface measurements is in principle the 
main aim of this thesis. 



2.2 Methods of measurements of the large-scale velocity fields 

Basically, there are three methods of measuring the photospheric velocity flelds: 
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Direct Doppler measurement - provides only one component (line-of-sight) of the 
velocity vector. These velocities are generated by local photospheric structures, am- 
plitudes of which are significantly greater than amplitudes of the large-scale velocities. 
The complex topology of such structures complicates an utihsation for our purpose. 
Analysing this component in different parts of the solar disc led to very important 
discoveries (e.g. supergranulation - Hart 1956 and Leighton et al. 1962). 



Tracer-type measurement - provides two components of the velocity vector. When 
tracing some photospheric tracers, we can compute the local horizontal velocity vec- 
tors in the solar photosphere. Tracking motions of sunspots across the solar disc led 



to the discovery of the differential rotation (Carrington, 1859). 

Local helioseismology - provides a full velocity vector. The local helioseismology (see 



Kosovichev 


1996 


Zhao et al. 


2001 



or 



Zhao & Kosovichev , 2004 ) is a very promising 



method using the information about the solar oscillations to infer the structure and 
also the dynamics in the convection zone. 

In this work we used mostly the tracer type method, local correlation tracking in particular, 
applied to the surface structures. 

The method needs a tracer - a significant structure recorded in different frames, the 
lifetime of which is much longer than the time lag between the correlated frames. We 
decided to use the supergranulation pattern in the full-disc dopplergrams, acquired by the 



Michelson Doppler Imager (MDI; Scherrer et al. , 1995) onboard Solar and Heliospheric 
Observatory (SoHO). We assume that supergranules are carried as objects by the large- 
scale velocity field. This velocity field is probably located beneath the photosphere, so that 
the resulting velocities will describe the dynamics in both the photospheric and subpho- 
tospheric layer. The existence of the supergranulation on almost the whole solar disc (in 
contrast to magnetic structures) and its large temporal stabihty make the supergranulation 
an excellent tracer. 



2.2.1 Local correlation tracking 

Since the photosphere is a very thin layer (0.04 % of the solar radius), the large-scale pho- 
tospheric velocity fields have to be almost horizontal. Then, the tracer-type measurement 
should be sufficient for mapping the behaviour of such velocities. In this field the local 
correlation tracking (LCT) method is very useful. 

This method was originally designed for the removal of the seeing-induced distortions 



in image sequences (November, 1986) and later used for mapping the motions of granules 



in the series of white-light images (November & Simon, 1988) under the name local cross- 



correlation. The method works on the principle of the best match of two frames that record 
the tracked structures at two different instants. 

The algorithm is applicable to two frames Ji and I2 having the same dimension that 
were captured in different instants. There is a time lag r between both frames, which 
has to be smaller than the lifetime of used tracer. For every pixel in the first frame a 
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Figure 1: Simple scheme explaining the performance of a LCT technique. 

subframe (correlation window) is chosen, described by the coordinates of its centre (xq, Uq) 
and the size p. Since in this study we used the Gaussian-weighted window, p is equal to 
the FWHM of the Gaussian profile. The parameter p is selected according to used tracer, 
so that p is larger than a characteristic size of the tracer. Let the subframe in image Ji is 

Si{xq, Uq) is compared with a subframe in image I2 of the same size, which has its centre 
in the coordinates (xq + 5x, Uq + 5y). Let the subframe in image I2 is 5*2 (xq + 5x, Uq + 5y). 
The proper motion of tracers in the point (xo,?/o) is defined by a displacement {5x,5y), 
which maximise the correlation of Si and 82- See Fig. [1} 

In this work we defined the correlation function C{xq + 5x, i/q + 6y) as weighted absolute 
difference of Si and 5*2: 

C{xo,yo,Sx,5y) = W{p) ■ \Si{xo,yo) - S'2(xo + 5x,yo + 5y)\ , (4) 

where W{p) is the weighting function, which depends on the size of the correlation window. 
We chosed W{p) with a Gaussian shape given by 

W{x,y) = e-'^, (5) 

where a = 0.424661p and the formula is valid for x G (— 2p, 2p) and y G {—2p, 2p). 

The correlation of each pair of subframes is performed for various values of 6x, 6y. In 
this work we use three integer values 6x and 6y given by 

6x = 0, ±shift; 6y = 0, ±shift, (6) 

where shift is a free parameter of LCT program. 

In our case for every pixel of the image with coordinates (xo,|/o), where xo G {0,Nx) 
and yo G (0, Ny) {N^ and Ny describe dimensions of input image) we obtain a 3x3-matrix 
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M(xo,?/o) containing the values of correlation of subframes 5*1 and 5*2 when displaced by 
5x, 5y. 



M{xo,yo) 



Cixo,yo,-s,+s) 
Cixo,yo, -s,0) 
C{xo,yo, -s, -s) 



C(xo,?/o,0, +s) C{xo,yo,+s,+s) 
C{xo,yo,0,0) C(xo,?/o,+s,0) 
C{xo,yo,0,-s) C{xo,yo,+s,-s) 



(7) 



where s = shift. 

Extremal value of M (maximum in case of the correlation is used, minimum if the 
absolute difference of subframes is used) is equivalent to the best match between Si and 
82- Corresponding values of 6x and Sy mean the displacement of the tracer in point (xq, yo) 
of the input frames. To obtain the subpixel precision, a surface is fitted to the values of 



M and the extremum is found on the fitted surface. Some numerical tests (November & 



Simon , 1988 ) showed that the best results are obtained by a biquadratic surface given by 

f{x, y) = ai + a2X + a^y + a^x'^ + a^y"^. (8) 



According to November & Simon ( 1988 ) the polynomial surface of order less than 2 over- 



estimate the displacements, while polynomial surface of order larger that 2 underestimates 



the displacements significantly. After some algebra (Darvann, 1991) we obtain for coeffi- 
cients a of ([s]) fitted to the matrix ([T]) formulae 

a, = ^M[0,0] + ^M[2,0] + ^M[0,2] + ^M[2,2] + ^M[0,1] + ^M[2,1]- 

-M[l,0]-M[l,l]-M[l,2], 
a2 = ^M[0,0] + ^M[2,0] + ^M[0,2] + ^M[2,2] + ^M[1,0] + ^M[1,2]- 

-M[0,l]-A/[1,1]-M[2,l], 

03 = M[0,0] -M[2,0] -M[0,2] +M[2,2], 

04 = -M[0,0 + M[2,2] + M[2,0] -M[0,1] +M[2, 1] -M[0,2], 
as = -M[0,0 + M[2,2] + M[1,0] -M[2,0] +M[0,2] -M[l,2]. 

Extremal shift of the best fit can be calculated using 



(9) 



min 



y 



1 
1 

15 



0-30,5 ~^ 7^02a4 



(10) 



where -D = | (yaia2 — a^a^j. The formulae assume that the absolute value of difference 
of subframes is used as the measure of correlation. The displacement of the tracers in point 
{xo,yo) is then given by 



= mirix ■ shift 
dy = miuy ■ shift. 
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The extrapolation in matrix M makes the resulting displacements unreliable, so that the 
displacements have to satisfy the conditions G {—shift, +shift), dy G {—shift, +shift). 
The desired behaviour can be accomplished by a suitable choice of free parameters of the 
LCT program. 

The result of LCT application to two frames is a two-dimensional map containing in 
each pixel two components of the displacement vector d = {dx,dy). The vector velocity 
field is obtained using 

(12) 



r 



For the evaluation of LCT in a time series of dopplergrams we used the modified 
program f lowmaker .pro implemented in IDL by Molowny-Horas & Yi (1994). 

LCT was recently used for tracking many features in various types of observations, 
especially for tracking the granules in high-resolution white-light images 



le.g. 



Sobotka 



et al. , 1999, Sobotka et al. , 2000). The same method was recently used for calculation 



Ambroz 


2001a 


and 


Ambroz 



2001b). Chae et al. (2004) measured the hehcity injected to a flaring active region using 



LCT applied to MDI full-disc magnetograms. Welsch et al. (2004) modified the LCT algo- 
rithm introducing the induction equation, which allows to determine all three components 
of surface flow field when tracking photospheric magnetograms. Smith et al. (2003) used 



LCT method for measurement of the chromospheric rotation profile when tracking full-disc 
Ha observations obtained at the Big Bear Solar Observatory. 

In some studies it has been shown that LCT technique underestimates the real veloc- 
ities due to the smoothing of processed data by a correlation window. For example in 



Georgobiani et al. (2007) the underestimation by 33 % is found. Georgobiani's study was 



done using different LCT code applied to simulated data and the results were compared 
with time-distance method applied to /-mode. This means that the correction factor is 
different for different settings of LCT method, and, therefore, should be determined (cal- 
ibrated) empirically for each particular study. Application of LCT to MDI Dopplergrams 
by jDeRosa fc Toomre (2004) showed the underestimation by more than 30 %. Some other 



studies also denote the underrepresenting of magnitudes by LCT (jSobotka et al. 
20 



Roudier et al. 1999 - 25 



1999 
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3 Method and tests on synthetic data 



A recent experience with applying this method to observed data (e. g. Svanda et al. , 2005 ) 
has shown that for the proper setting of the parameters and for the tuning of the method, 
synthetic (model) data with known properties are needed. The synthetic data for the analy- 
sis come from a simple numerical simulation (SISOID code = Simulated ^upergranulation 
as Observed In Dopplergrams), with the help of which we can reproduce the supergra- 
nulation pattern in full-disc dopplergrams. The design of the synthetic dopplergrams with 
known parameters is very important for the calculation of the vector velocity fields, so that 
we carried out the simulation carefully to get the correct and valuable results about the 
abilities of the method. 



3.1 The SISOID code 



The SISOID code is not based on physical principles taking place in the origin and evolution 
of supergranulation, but instead on a reproduction of known parameters that describe 
the supergranulation. Individual synthetic supergranules are characterised as centrally 
symmetric features described by their position, lifetime (randomly selected according to the 
measured distribution function of the supergranular lifetime - DeRosa et al. 2000 see Fig. [2] 
left), maximal diameter (randomly according to its distribution function that is basically 
normal with the mean of 31 200 km and the variation of 2 300 km - Wang & Zirin , 1989[ ) and 
characteristic values of their internal horizontal and vertical velocity components (randomly 
according to their distribution function with a normal shape described by parametres 
< Vh >= 258 ± 1 ms-\ <v^>=29±2 ms"^ - Hathaway et al. , 2002). The magnitude 
of the horizontal and vertical component of the internal velocity field depends only on the 
relative distance from the centre of the cell and is aproximated by curves shown in Fig. [3] 
The most important simplification in the SISOID code is that individual supergranules 
do not influence each other, but simply overlap. The final line-of-sight velocity at a certain 
point is given by the sum of line-of-sight velocities of individual synthetic supergranules at 
the same position. 

New supergranules can arise inside the triangle of neighbouring supergranules (iden- 



tification of such triangles is done by the Delaunay triangulation algorithmed by Barry 



1991) only when the triangle is not fully covered by other supergranules and when any of 



the supergranules located at the vertices of the triangle is not too young, so that in the 
future it could fully cover the triangle. The position of the origin of the new supergranule 
is the centroid of the triangle; each vertex is weighted by the size of its supergranule. The 
diameter evolves according to Fig. |2] right - in the lifetime the supergranule grows from 
zero (during first 40 per cent of its lifetime) to its maximum diameter and shrinks to zero 
again (during last 60 per cent of its lifetime). All the supergranules behave in this way, 
which roughly approximates the real behaviour of the convection cells. 

■k This chapter was in a condensed form pubHshed in Svanda, M., Klvaha, M., and Sobotka, M., 2006, 
Large-scale horizontal flows in the solar photosphere. I. Method and tests on synthetic data, Astronomy 
and Astrophysics, 458, 301-306. 
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Figure 2: Left - The histogram of supergranule lifetime used in the SISOID code. Adopted from 



DeRosa et al. (2000). Right - Approximated evolution of the diameter of synthetic supergranule 



during its lifetime as used in the SISOID code. 
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Figure 3: The dependence of the values of horizontal ('left ) and vertical ( right ) components of 
the internal supergranular velocity field on the radial distance from the centre of the synthetic 
supergranule, used in the SISOID code. 



The SISOID simulation is done in the pseudocylindrical Sanson-Flamsteed coordinate 



grid (Calabretta & Greisen, 2002); the transformation from the heliographic coordinates 



is given by the formulae: 

x = 'dcos({), y = (p, (13) 

where x and y are coordinates in the Sanson-Flamsteed coordinate system, and ^ and if are 
heliographic coordinates originating at the centre of the disc. At each step an appropriate 
part of the simulated supergranular field is transformed into heliographic coordinates. 
The output of the program is a synthetic dopplergram of the solar hemisphere in the 
orthographical projection to the disc. We assume in our simulation that the Sun lies in 
an infinite distance from the observer and that the P (position angle of the solar rotation 
axis) and bo (heliograhic latitude of the centre of solar disc) angles are known. 
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Each step in calculation includes the evaluation of the parameters of individual super- 
granules, then small old supergranules under the threshold (2 Mm in size) are removed 
from the simulation and all the triangles are checked, whether a new supergranule can 
arise inside them (see Fig. |4]). This step in the SISOID code corresponds to 5 minutes in 
real solar time. The computation is always started from the regular grid. Properties of 
"supergranules" are chosen randomly according to their real distribution functions. The 
first 1000 steps are "dummy", i. e., no vector velocity field is included and no synthetic 
dopplergram is calculated. This starting interval is taken for the stabilisation of the super- 
granular pattern. In the next steps, the model vector velocity field is already introduced. 
This field influences only the positions of individual cells. The dopplergram is calculated 
every third step. For one day in real solar time, 96 dopplergrams are calculated. 

The model velocity field with Carrington rotation added is applied according to the 
assumption of the velocity analysis that the supergranules are carried by a velocity field on 
a larger scale. In the simulation, the position of individual supergranules is influenced, and 
no other phenomena are taken into account. These synthetic dopplergrams are visually 
similar to the real observed dopplergrams (see Fig. [s]). 



3.2 Method of data processing 

The MDI onboard SoHO acquired the full-disc dopplergrams at a high cadence in certain 
periods of its operation - one observation per minute. These campaigns were originally 
designed for studying the high-frequency oscillations. The primary data contain lots of dis- 
turbing effects that have to be removed before ongoing processing: the rotation line-of-sight 
profile, p-modes of solar oscillations. We detected some instrumental effects connected to 
the data-tranfer errors. It is also known that the calibration of the MDI dopplergrams 
is not optimal (e.g. Hathaway et al. , 2002) and has to be corrected to avoid systematic 



Start from a regular grid 



After 1 step 



After 1000 steps 



Final synthetic dopplergram 



Properties ol individual supergranules are chosen 
rsnOomly according to ifieir measured distribution function 



First 1000 steps are 'dummy" to 
stabilise the supergranular pattern 




Points where the new 
synthetic supergranule wilt be 
formsd In the next step 



Neightjouring supergranules do 
not influence each other, they 
simply overlap; final lirte-ot-sight 
velocity is given by the sum of 
individusl line-of-sight velocities 




Figure 4- The cartoon describing how the SISOID code works and forms the final synthetic 
dopplergram 
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Figure 5: Comparison of the real dopplergram ('left, observed by MDI/SoHO) and the simulated 
one ('rightj. Both images are visually very similar. The black colour means line-of-sight velocity 
—700 ms~^ (towards an observer) while white represents +700 ms^^. 



errors. While examining long-term series of MDI dopplergrams, we have met systematic 
errors connected to the retuning of the interferometer. We should also take those geomet- 
rical effects into account (finite observing distance of the Sun, etc.) causing bias in velocity 
determination. According to Strous (2000) for example, the bias coming out of a perspec- 
tive is about 2 ms~^, and it depends on the position on the disc. It has been proven (e. g. 



Liu & Norton, 2001) that MDI provides reliable velocity measurements when the magnetic 
field is lower than 2000 Gauss. The velocity observation by MDI will induce up to 100 % 
error if the magnetic field is higher than 3000 Gauss due to the magnetic sensitivity of 
the used Nil line and the limitations of computational algorithm, which cause crosstalks 
between measured MDI dopplergrams and magnetograms. The removal of these effects 
will be described in detail in Section |5| while the synthetic data used in this study do not 
suffer from these phenomena. 

As input to the data processing we take a one-day observation that contains 96 full-disc 
dopplergrams in 15-minute sampling. Structures in these dopplergrams are shifted with 
respect to each other by the rotation of the Sun and by the velocity field under study. The 
primary data must be pre-processed by removal of the manifestation of the Carrington 
rotation and by the suppression of the p-modes. 

First, the shift caused by the rotation has to be removed. For this reason, the whole data 
series (96 frames) is "derotated" using Carrington rotation rate, so that the heliographic 
longitude of the central meridian is equal in all frames and also equals the heliographic 
longitude of the central meridian of the central frame of the series. This data-processing 
step causes the central disc area ("blind spot" caused by prevailing horizontal velocity 
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Figure 6: The scheme showing individual steps during the preparation of the dopplergrams for the 
LCT technique. 



component in super granules) in the derotated series to move with the Carrington rate. 
During the "derotation" the seasonal tilt of the rotation axis towards the observer (given 
by 6o ^ heliographic latitude of the centre of the disc) is also removed, so that 6o = in all 
frames. 

Then the data series is transformed to the Sanson-Flamsteed coordinate system to 
remove the geometrical distortions caused by the projection of the sphere to the disc. 
Parallels in the Sanson-Flamsteed pseudocylindrical coordinate system are equispaced and 
projected at their true length, which makes it an equal area projection. Formulae of the 



transformation from heliographic coordinates are given by Eq. (13). 

The noise coming from the evolutionary changes in the shape of individual supergranules 
and the motion of the "blind spot" in the data series with the Carrington rotational rate 



are suppressed by the k-uj filter in the Fourier domain (Title et al. , 1989, Hirzberger et al. 
1997). The cut-off velocity is set to 1500 ms^^ and has been chosen on the basis of 
empirical experience. The procedure can be seen in Fig. [6j 

The existence of the differential rotation complicates the tracking of the large-scale 
velocity field, because the amplitudes and directions of velocities of the processed velocity 
field have a significant dispersion. We have found that, when the scatter of magnitudes 
is too large, velocities of several hundred ms~^ cannot be measured precisely by the LCT 
algorithm where the displacement limit for correlation was set to detect velocities of several 
tens of ms~^. Therefore the final velocities are computed in two steps. The first step 
provides a rough information about the average zonal flows using the differential rotation 
curve 

uj = A + Bsm^h + CsmH , (14) 

and calculating its coefficients. 

In the second step this average zonal flow is removed from the data series, so that during 
the "derotation" of the whole series the differential rotation inferred in the first step and 



expressed by ( 14 ) is used instead of the Carrington rotation. The scatter of the magnitudes 
of the motions of supergranules in the data transformed this way is much smaller, and a 
more sensitive and precise tracking procedure can be used. 
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The LCT method is used in both steps. In the first step, the checked range of velocity 
magnitudes is set to 200 ms~\ but the accuracy of the calculated velocities is roughly 
40 ms~^. In the second step the range is only 100 ms~^ with much better accuracy. The 
lag between correlated frames equals in both cases 16 frame intervals (i. e. 4 hours in solar 
time), and the correlation window with FWHM 30 pixels equals 60" on the solar disc in 
the linear scale. In one observational day, 80 pairs of velocity maps are calculated and 
averaged. 

For the calculation we use the adapted program f lowmaker .pro originally written in 
IDL by Molowny-Horas & Yi (1994). The algorithm has a limitation in the range of 
displacements that are checked for each pixel. The quality of correspondence (in our case 
the sum of absolute differences of both correlation windows) is computed in nine discrete 
points, then the biquadratic surface is fitted through these nine points, and an extremum 
position (Darvann, 1991) is calculated (see Section 2.2.1). The final displacement vector 
is equivalent to the position of the extremum. 



3.3 Results of the synthetic data experiment 

In our tests we have used lots of variations of simple axisymmetric model flows (with 
a wide range of values of parameters describing the differential rotation and meridional 
circulation) with good success in reproducing the models. When comparing the resulting 
vectors of motions with the model ones, we found a systematic offset in the zonal component 
equal to foSset, zonal = "15 ms^^. This constant offset appeared in all the tested model 
velocity fields and comes from the numerical errors during the "derotation" of the whole 
time series. For the final testing, we used one of the velocity fields obtained in our previous 



work (Svanda et al. , 2005). This field approximates the velocity distribution that we may 
expect to observe on the Sun. The model flows have structures with a typical size of 60", 
since they were obtained with the correlation window of this size. 

The calculated velocities (with f offset, zonal = —15 ms~^ corrected) were compared with 
the model velocities (Fig. [?]). Already from the visual impression it becomes clear that 
most of vectors are reproduced very well in the direction, but the magnitudes of the vec- 
tors are not reproduced so well. Moreover, it seems that the magnitudes of vectors are 
underestimated. This observation is confirmed when plotting the magnitudes of the model 
vectors versus the magnitudes of the calculated vectors (Fig. [sjleft). The scatter plot con- 
tains more than 1 million points, and most of the points concentrate along a strong linear 
dependence, which is clearly visible. This dependence can be fitted by a straight line that 
can be used to derive the calibration curve of the magnitude of calculated velocity vectors. 
The calibration curve is given by the formula 

Vcor = 1.13fcalc, (15) 

where fcaic is the magnitude of velocities coming from the LCT, and fcor the corrected 
magnitude. The directions of the vectors before and after the correction are the same. The 
uncertainty of the fit can be described by 1-cr-error 15 ms~^ for the velocity magnitudes 
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Figure 7: Left - model vector velocity field. Right - a velocity field that was computed by applying 
the LCT m,ethod to the synthetic supergranulation pattern with the imposed model field. The arrow 
lengths, representing the velocity magnitudes, have the same scale. The images are visually very 
similar, however the magnitudes of calculated velocities are underestimated. 
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Figure 8: Left - Scatter plot for inference of the calibration curve. Magnitudes of calculated 
velocities are slightly underestimated by LCT, but the linear behaviour is clearly visible. A line 
representing the 1:1 ratio is displayed. The calibration affects only the magnitudes of the flows, 
while the directions do not need any correction. Right - Dependence of the 1-a-error of the 
calculated velocity on the time lag when no model velocity field was introduced and only the 
evolution of supergranules have been taken into account. The dashed line represents lag 16 (4 
hours) that is usually used in our method. 
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Figure 9: Influence of the calculated velocity field on the choice of FWHM of the correlation 
window. Left - 12Cf' , right - 200" . The model velocity field is the same as in Fig. 
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Figure 10: Histograms of velocity magnitudes for various FWHM of the LCT algorithm. 



under 100 ms~^ and 25 ms~^ for velocity magnitudes greater than 100 ms^^. The uncer- 
tainties of approx. 15 ms~^ have their main origin in the evolution of supergranules. We 
studied the dependence of the error of velocity determination on the time lag used when 
no model velocity field was introduced. We found that this dependence is slowly increasing 
with the time lag (Fig. [s] right) due to the evolution of individual supergranules. Evolution 
of supergranules is only one part of story, but it gives the lower limit of accuracy that can 
be obtained by this method. We also ran a test of the LCT sensitivity on the evolution 
of supergranules when a known underlying velocity field is introduced and came to similar 
results. 
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We tested the sensitivity of the method to the choice of values of FWHM of the cor- 
relation window and of the lag between correlated frames in the LCT method. We found 
that our method is practically insensitive to the choice of the time lag between correlated 
frames when the lag is in the interval of 10-24 (2.5-6 hours). The larger the lag we choose, 
the lower the velocities we are able to detect. On the other hand, we have to take into 
account that a larger lag between correlated frames causes more noise in calculated results 
coming from evolutionary changes of supergranules and probably also from evolutionary 
changes in the velocity field under study. According to our tests, for a time lag greater 
than 30 (7.5 hours), the numerical noise raises very fast. The lag 16 (4 hours) seems to be 
a good tradeoff between sensitivity and noise. 

The choice of different FWHMs of the correlation window changes the spatial resolution 
according to FWHM. The general character of the vector field is preserved within the 
limits of resolution (cf. Fig. |9]). Larger FWHM causes a smoothing of results and an 



underestimation of vector magnitudes (cf. Fig. 10 ). We found the used parameters (FWHM 
60", lag 4 hours) to be the best compromise, however these values can be changed during 
the work on real data. 

Using the synthetic data generated by the SISOID code we have verified that the 
proposed method is reliable when measuring the large-scale velocity fields in the solar 
photosphere using the LCT method applied to full-disc dopplergrams. 
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4 Direct comparison to the time-distance helioseismology 



Velocity fields in one instant may be calculated using different techniques. However, the 
results obtained by different methods may have discrepancies. These can be caused by the 
nature of the methods, e. g. due to different types of averaging, and also because of the 
use of different datasets from different instruments. In addition, various disturbing effects 
can be important. Therefore, we decided to compare the results obtained by two different 
methods, time-distance helioseismology and local correlation tracking (LCT), using the 
same set of data: high-cadence dopplergrams covering almost one Carrington rotation 



obtained from Michelson Doppler Imager (MDI, Scherrer et al. 1995) 



Both methods provide surface or near-surface velocity vector fields. However, the results 
of these methods can be interpreted differently. While local helioseismology measures 
intrinsic plasma motions (through advection of acoustic waves), LCT measures apparent 
motions of structures (granules or magnetic elements). It is known that some structures do 
not necessarily follow the flows of the plasma on the surface. For example, supergranulation 



appears to rotate faster than the plasma (Beck & Schou, 2000), which may be caused 



by travelling waves (Gizon et al. 2003) or may be explained also as projection effect 



(Hathaway et al. , 2006 ). Some older studies (see e. g. Rhodes et al. , 1991 ) also suggest that 



the difference in flow properties measured on the basis of structures motions and plasma 
motions is caused by deeper anchor depth of these structures. An evolution of pattern 
may also play significant role (e. g. due to emergence of magnetic elements). Another 
possibility is that surface structures are not coherent features, but patterns traveling with 
a different group velocity than the surface plasma velocity, such as occurs for the features 



present in simulations of travelling- wave convection (e. g. Hurlburt et al. , 1996). 

Some attempts to compare the results of local helioseismology and the LCT method for 
large scales, with characteristic size 100 Mm and more, have been carried out by Ambroz 



(2005), but his results were inconclusive. The correlation coefficient describing the match 



of the velocity maps obtained by local helioseismology and the LCT method was close 
to zero. Nevertheless, there were compact and continuous regions of the characteristic 
size from 30 to 60 heliographic degrees with a good agreement between the two methods, 
so that one could not conclude that the results were completely different. In his study 
many factors could play a significant role: the techniques were applied to different types 
of datasets (LCT was applied to low resolution magnetograms acquired at Wilcox Solar 
Observatory, and the time-distance method used MDI dopplergrams). Both techniques 
had very different spatial resolution, and also the accuracy of the measurements was not 
well known. 

We decided to avoid these problems and analyze the same data set from the MDI 
instrument on SoHO. The MDI provides approximately two months of continuous high- 
cadence (1 minute cadence) full-disc dopplergrams each year. This Dynamics Program 
provides data suitable for helioseismic studies, and also for the local correlation tracking of 

★ This chapter was in the condensed form pubhshed as Svanda, M., Zhao, J., and Kosovichev, A. G., 
2007, Comparison of Large-Scale Flows on the Sun Measured by Time-Distance Helioseismology and Local 
Correlation Tracking, Solar Physics, 241(1), 27-37. 
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supergranules. Thus, this is a perfect opportunity to compare performance and results of 
two different techniques using the same set of data, and avoid effects of observations with 
different instruments or in different conditions. 



4.1 Data preparation 

The selected dataset consists of 27 data-cubes from March 12th, 2001, 0:00 UT to April 6th, 
2001, 0:31 UT, where each third day was used, and in these days three 8.5-hour long data- 
cubes were processed (so that every third day in the described interval was fully covered by 
measurements). Each data-cube is composed of 512 dopplergrams (with spatial resolution 
of 1.98" px~^) at a one-minute cadence (so that covering 8 hours and 32 minutes). All the 
frames of each data-cube were tracked with a rigid rate of 2.871 yurads"^, remapped to the 
Postel projection with a resolution of 0.12 ° px~^ (1,500 kmpx~^ at the center of the disc), 
and only a central meridian region was selected for the ongoing processing (with size of 
256x924 px covering 30 heliographic degrees in longitude and running from —54° to -|-54° 
in latitude), so that effects of distortions due to the projection do not play a significant 
role. 

Tracked data-cubes were used to perform the time-distance analysis. From all the 
frames in each data-cube the mean dopplergram (like Fig. [iSjleft) was subtracted to sup- 
press the influence of velocity structures like supergranulation and to highlight the signals 
of p-modes of solar oscillations. The surface gravity wave (/-mode) has different disper- 
sion characteristics than p-modes used in this study, and, therefore, it is filtered out from 
the k-uj diagram before computing the travel times. The /-mode, if not filtered out, will 
disturb p-modes measurements, and it is also not straightforward to perfrom inversions if 
not separating two different modes. The p-mode inversions are less sensitive to the surface 
flows than /-mode data, but still recover the large-scale flows well. 

P-modes of solar oscillations have their origin in the solar convection zone and travel 
through the solar interior to the surface. The time of the excursion of the wave bulk 
depends on the speed of sound and on the velocity field describing the mass flow in the 
layers of the solar interior, through which the bulk is travelling. In the time-distance 
technique, the travel times of oscillations bulk from the point in the photosphere (central 
point) to surrounding annuli around this point are measured. The radius of the studied 
diameter of the annuli is related to the depth in the interior, where the studied oscillation 
mode is reflected back to the surface. Travel times are measured by the cross-correlation 
between Doppler velocities in the central point and velocities in the selected annuli around 
this point. 

The mass flow velocities in the interior are calculated from the differences of travel times 
from the central point to the surrounding annuli and the travel times from the surrounding 
annuli to the central point when the state properties in the affected layers of the solar 
interior are known. In this study the theoretical travel ray approximation is derived from 



the solar model S ( Christensen-Dalsgaard et al. , 1996). Dividing the annuli into sectors the 



underlying flow field of selected orientations can be inferred. For details see Kosovichev 



(1996 


), 


Zhao et al. 


(2001 


), or 
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The time-distance inversion results were smoothed by a Gaussian with FWHM of 30 px 
to match the resolution to the LOT method, and only the horizontal components {v^, Vy) 
of the full velocity vector were used. 

While for the time-distance method the p-modes of solar oscillations play a crucial role, 
they significantly influence the performance of the LCT method in a negative way. The 
oscillations are clearly visible in the dopplergrams, and, thus, cause random errors in the 
calculation of displacements. Therefore, before applying the LCT method the oscillation 
signals must be suppressed. For our high-cadence data it is possible to do this using 



temporal averaging. According to Hathaway (1988) it is better to use a Gaussian type of 



temporal averaging than the boxcar one. We average the dopplergrams over 31 minute 
periods with weights given by a formula 
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where At is the time between a given frame and the central one (in minutes), 6 = 16 min 
and a = 8 min. We verified that this filter suppresses the solar oscillations in the 2-4 mHz 
frequency band by a factor of more than five hundred. 

The other issue significantly influencing the performance of the LCT method is the 
change of contrast and background intensity caused by solar rotation. Due to tracking the 
Doppler images with rotation, the magnitude of the line-of-sight component of the solar 
rotational velocity changes from frame to frame and affects the LCT results. The method 
interprets these changes as motion towards east, mainly in the central part of the solar 
disc, where the contrast in the structures of dopplergrams is very low (see Fig. |15| left). We 
suppress the influence of the moving background by subtraction of a polynomial surface fit 
of the third order. We have tested that this provides almost the same results as the other 
possible procedures: local removal of the mean values and unsharp masking. Subtraction 
of the polynomial fit is not so sensitive to anomalies in the dopplergrams, caused by regions 
with strong magnetic field. 

The LCT method used in this study is described by the following parameters: the 
time-lag between correlated frames is 120 frames (2 hours), the correlation window has a 
Gaussian shape with FWHM of 30 px, the correlation is measured by the sum of absolute 
differences of subframes (it is faster than calculation of the correlation coefficient and 
provides the same results), the extremum position is calculated using the nine-point method 



of Darvann (1991). For each data-cube, the results of all the correlated pairs are averaged, 
so that the method provides an averaged flow field in 8.5 hours in the same sense as the 
time-distance analysis. 



4.2 Results 

4.2.1 Statistical processing 

The results containing 27 horizontal flow fields were statistically processed to obtain the 



cross-calibration curves for these methods. It is generally known (see discussion on page 13) 
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Figure 11: Left -Vx^ lot versus v^, t-d plot. Red crosses denote mean zonal velocities (differential 
rotation) which have been used for fitting of the calibration formula |T7[ j. Right - the same for 
Vy component, the regression fit is described by Eq. (18). 



that the LCT method sUghtly underestimates the velocities; thus, the results should be 
corrected by a certain factor. From the comparison of the x-component of velocity (cf. 
Fig [ll] left) we obtained parameters of a linear fit given by (numbers in parentheses denote 
a l(T-error of the regression coefficient) 

Vx, LCT = 0.895(0.008)^;,, t-d - 12.6(0.3) m s'^ (17) 

The correlation coefficient between Vx, lct and Vx, t-d is p = 0.80. We assume that the 
time-distance measurements for Vx, t-d are correct, and the magnitude of the LCT measure- 



ments, Vx, LCT, niust be corrected according to the slope of Eq. (17). This correction factor 
has a value of 1.12, which is in perfect agreement with the correction factor of 1.13 found 
in the tests of the same LCT code using synthetic dopplergrams with the same resolution 
and similar LCT parameters (see Section [3]). We assume that both velocity components 
obtained with the LCT method should be corrected by this factor. 



The regression line of Vy component (Fig. 11 right) is 

Vy, LCT = 0.56(0.01)^;^,, t-d + 0.4(0.2) ms^^ 
After the slope correction using the Vx fits, the regression curve is slightly different: 



V. 



, LCT = 0.63(0.01)t;j,, t-d + 0.4(0.2) ms-\ (19) 



with the correlation coefficient between Vy, lct and Vy, t-d close to 0.47. The slope of the 
linear fit differs significantly from the expected value 1.0. For this study we decided to 
correct the y-component of the time-distance results. 

We have tested that this asymmetry is not related to the LCT technique. The tests did 
not show any preference in direction of flows measured by LCT or any dependence of the 
results on the size of the fleld of view (which does not have a square shape in our case). 
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Figure 12: Histogram of the angular differences (^if) between directions of the velocity vectors 
obtained by the time-distance and LCT techniques. Dashed vertical line denotes the mean value 
and solid vertical line represents the median of A99. 



Also in the study performed in Section |3] based on synthetic data the asymmetry between 
the zonal and the meridional component was not encountered. 

We found two possibilities that could explain this behaviour. The first explanation 
is a drift of the supergranular pattern towards the equator. This case does not explain 
why the meridional velocities from both techniques seem to be proportional to each other. 
A systematic drift would rather be depicted as a systematic constant shift, or a shift 
depending on the latitude. However, the meridional components of velocities are generally 
rather small, so that the errors of the measurements can play a significant role and the 
proportional behavior can be only apparent. 

The second explanation is based on unspecified asymmetries influencing travel-time 
measurements, for instance, due to different sensitivity of the MDI instrument to p-modes 
propagating in the east-west and north-south directions. As it has been studied recently 



(Georgobiani et al. , 2007), a comparison between /-mode time-distance and LCT applied 
to realistic numerical simulation did not show such asymmetry. The asymmetry in the 
east-west and north-south directions observed by time-distance helioseismology was on the 



contrary noticed in a recent study based on numerical simulated data (Zhao et al. , 2007). 



So that the asymmetry takes place only in the p-modes inversions and should be further 
investigated in more details. 

The final calibration formulae providing the best statistical agreement between the 
velocities calculated using both methods are: 
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t'z.LCTxorr — 1- 12f 2:,LCT,calc (20) 

,corr (21) 

^x,t— d,corr ^x,t— d,calc (22) 

0.63fy,t-d,calc, (23) 

where the index corr denotes the corrected value, and the index calc denotes the original 
calculated value. 



After the corrections, as presented in the histogram of Fig. [12} the differences between 
the directions of the velocity vectors (Aip) calculated by these techniques are quite reason- 
able. The mean value of the distribution is 43.56°, however, the mean value is not a good 
indicator in this case because the distribution function is not normal. The median value is 
24.02°; and 66.6 % points have the difference in the corresponding vector directions under 
45°. 

Instead of computation of the correlation coefficient of the arguments of both vector 
fields we decided to compute a magnitude-weighted cosine of Aip. This quantity is given 
by 



\vt-d-VLCT\ 



PW = , , (24 j 

\^t-d\ 

where v^.^ is the time-distance vector field, "Vlct is the LCT vector field and the summation 
is performed over all vectors in the field. The closer this quantity is to 1, the better is the 
agreement between both vector fields. Larger vectors are weighted more than smaller ones. 
We have found that in our case pw = 0.86, which means an almost perfect match. 

4.2.2 Mean velocities 

In addition to the detailed comparison of the vector fields, we compare the mean flows, 
the differential rotation and the meridional circulation. These flows can be quite simply 
calculated from the results of both techniques. In both cases, they provide the mean zonal 
and mean meridional flows for the Carrington rotation No. 1974. The results are displayed 



in Fig. 13 , where the differential rotation curves are compared with a standard profile from 



Snodgrass & Ulrich (1990) in the left panel. It can be clearly seen that the latitudinal 
profiles for both techniques are very similar and also that the mean velocities do not differ 
so much in magnitude. The correlation coefficients are p = 0.98 for the zonal flow and 
p = 0.88 for the meridional flow. In the differential rotation curves, the LCT results give 



a little slower rotation, which is also seen from Eq. (17). The mean difference of average 
zonal velocities obtained by both techniques is 14.1 ms~^ (see Fig. 14). 

It is coherent with the systematic shift found in Section [3} i. e. a systematic offset of 
—15 ms^^ caused by the change of the projection effect due to the "derotation" of the 
studied data series. We believe that the offset found between the LCT and time-distance 
results is caused by the same reason. 
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Differential rotation 




Figure 13: Left: Thick curve shows the mean zonal velocity as a function of latitude (differential 
rotaton) obtained from the time-distance data, thin curve shows the mean zonal velocity obtained 
from the LCT data. Dashed line represents a standard rotation profile ( "Snodgrass rate"). The 
mean zonal velocities are plotted in the coordinate system rotating rigidly with 2.871 firads^^. 
Right: Velocities of the mean meridional flow as a function of latitude obtained by the time- 
distance technique (thick curve) and by the LCT method (thin curve). The correlation coefficient 
for the mean zonal flow is p = 0.98 and for the mean meridional flow p = 0.88. 




b [deg] 



Figure 14: Mean zonal flows (like in Fig. 13 left) with the standard rotation profile ("Snodgrass 



rate") removed. The systematic shift between the curves obtained from time-distance results (thick 
line) and LCT results (thin line) is evident and is in average I4.I ms^^ . 



We conclude that for the mean flows the results obtained by two such different tech- 
niques agree very well. 



4.2.3 Detailed comparison 

For a detailed comparison of the flow fields, we selected one data cube, representing 8.5- 
hour measurements centered at 4:16UT of March 24th, 2001 (/o=214.3°; see the averaged 
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Figure 15: Part of the solar disc chosen for the detailed comparison of the results of time- distance 
and LCT method. The 8.5-hour measurements are centered at March 24th, 2001, 4-16 UT. One 
pixel represents 0.12° in heliographic coordinates. Left - Averaged MDI dopplergram. Note the 
low contrast of the supergranular cells in the center of the images. This is the "blind spot" caused 
by the prevailing horizontal motions in the photosphere with this resolution. Right - Averaged 
MDI magnetogram. 



MDI Dopplergram and magnetogram in Fig. 15). In this map, the correlation coefficient 
for the a;-component of the velocity is p = 0.82 and for the y-component: p = 0.58, and 
for the vector magnitude: p = 0.73. The vector plots of the flow fields obtained by both 
techniques, shown in Fig. 16, seem to be quite similar, in general; however many differences 
can be seen. The regions where the differences are most significant correspond to relatively 
small (under 50 ms~^) velocities. This is clear from the map of the differences between 



the vector directions displayed in Fig. 17 



The statistics of the differences between the directions of the corresponding vectors 
(Ay^) is presented in Fig. 12 Values of Aip slightly anti-correlate with the averaged magni- 
tude of the corresponding vectors (p = —0.58). We think that this is due to uncertainties of 
both techniques. From our tests using synthetic data, it became clear that the inaccuracy 
of the local correlation tracking code is 15 ms~^ for velocities smaller than 100 ms~^ and 
25 ms~^ for velocities larger than 100 ms~^, for both components (see Section We 



think that the 10% accuracy for the time-distance velocity vectors is a reasonable esti- 
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Figure 16: Left - velocity field obtained by the LCT method. Right - vector field obtained by 
the time- distance technique. Both plots are centered at heliographic coordinates bo = 0.0°, Iq = 
214.3°, units on both axes are pixels in the data frame with resolution 0.12° px~^ in the PosteVs 
projection. 



mate. As Zhao et ah (2001) stated, cross-talk effects between horizontal flows and vertical 
flow components of flow velocities affect the time-distance inversion results. The cross-talk 
effect prevents us from inverting the vertical velocity correctly, but does not block the 
determination of horizontal velocities satisfactorily (Zhao & Kosovichev 2003). However, 
the vertical velocities are not discussed in this Section because they are not measured 
by the LCT technique. Obviously, the inaccuracy in one component may cause a signifl- 
cant change of the direction of the horizontal vector for small velocities; and, hence, the 
agreement of both techniques in such areas is not as good as in the areas of high velocities. 

The vector velocity fleld may be also influenced by the temporal evolution of the traced 
pattern. We have tested using the full-disc MDI dopplergrams that temporal changes at 
mesogranular and smaller scales are effectively flltered out by a k-uj fllter. The temporal 
evolution of the supergranular pattern may, in the worst case, signiflcantly influence the 
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Figure 17: Difference of phases of corresponding vectors (Aip). Left - for the full vector data, 
right - for the data with mean zonal flow subtracted. In the second case the magnitudes of vectors 
everywhere in the studied area are comparable and the inaccuracies of individual measurements 
equally important. With the mean zonal flow removed, the distribution of Aip becomes more 
uniform in the studied area. 



calculated vector velocity field in the close (roughly equal to the FWHM of the chosen 
correlation window) vicinity of a rapidly changing (e. g. disappearing) supergranule. 



4.3 Conclusions 

Flow velocity fields on the solar surface obtained by two different techniques, time-distance 
helioseismology and local correlation tracking (LCT), were compared. Despite the fact that 
the first technique uses p-modes of solar oscillations to compute the velocity field (in the 
data, the large-scale structures like supergranulation were suppressed), while the other one 
uses large-scale supergranulation pattern from averaged Doppler images as tracers for the 
velocity vectors determination (and requires p-modes removal), we found that both results 
match reasonably well. We have confirmed some recent studies that the LCT method 
slightly underestimates the actual velocities (as the consequence of a smoothing procedure), 
and we determined empirical correction factors. After the corrections, the match in global 
velocity structures, mean zonal and meridional flows, is very good. The results of a detailed 
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Figure 18: Left - Velocity field obtained by the LCT method applied to the data cube used in the 
helioseismic study. Right - Velocity field obtained by the 'standard' LCT method described in 
Section^^ One can see that velocity field obtained using the same LCT code and different data 
processing may differ. The magnitude weighted cosine has a value of = 0.71, so both fields are 
still similar. There are two factor in the data processing that cause the difference: Velocity field 
obtained from the helioseismic data cube are in the Postel projection and the velocity field using 
the 'standard' method was calculated using the lag of 4 hours and averaged over 24 hours, while 
velocity field obtained from the helioseismic data cube was calculated using the lag of 2 hours and 
averaged over 8.5 hours. 

comparison of the vector velocity fields are not so satisfactory. However, the correlation 
coefficient for individual components of the flows is positive and significant, so we conclude 
that a meaningful match is found. It is shown that the largest disagreement is caused 
by very small velocities in some regions, where the errors of both methods become quite 
significant. 
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5 Application to real full-disc data 

5.1 Processed data 

In agreement with the aim of processing of the full-disc MDI dopplergrams we have pro- 
cessed all the suitable dopplergrams measured in the Dynamic campaigns approximately 
two months each year since 1996. MDI provides in this way the most homogeneous mate- 
rial, which can be used for such studies. Similar data series are measured by the GONG 
network, but it is based on ground-based instruments and the measurements are strongly 
influenced by the turbulence in the Earth's atmosphere. 

Among all, the MDI data suffer of some issues. The gaps in the data series are the 
most important problem. In some years, the gaps were present almost every day and had 
a length of several hours. They are mostly caused by the misoperation of the instrument. 
Other important issues are caused by the errors during the transmission of the data from 
SoHO to the central data storage. They usually depict as missing parts of the frame. Some 
frames also have a wrong orientation. 

Since the homogeneous non-interrupted series are needed for the current study, we 
have to deal with the described errors. The missing parts of the frames and the issue of 
the wrong orientation are non-correctable, therefore we delete all the frames with missing 
values and inconsistent telemetry in headers of the frames. The gaps are filled with linearly 
interpolated frames to run the reduction routines smoothly. 

The linear interpolation of the missing frames influence the data processing in a neg- 
ative way. During the temporal averaging in order to suppress the p-modes of the solar 
oscillations, the series containing interpolated frames do not fulfill the assumptions of the 
processing. This in fact leads to the enhancement of the solar oscillations signature in 



dopplergrams - see Fig. 19 



Dopplergrams influenced by the linear interpolation are not used in the ongoing data 
processing routine, because the signal of p-modes confuses the LOT algorithm. The dif- 
ference between a "bad" dopplergram and a "good" one is seen mainly in the centre of 
the disc, where in the "good" dopplergram the signal of supergranulation is weak, while in 
the "bad" one the signal of p-modes is clearly visible on the centre of the disc. The "bad" 
dopplergram can be therefore easily detected on the basis of excess of power in the power 
spectra in the band of 0.32-0.34 px~^, which corresponds to the size of 3 px (4500 km) - 



see Fig. 20 



The processed data belonging to the Dynamics campaign are summarized in Tab. [Tj 
Each day, two 24-hours series separated by 12 hours were processed. If the series contains 
some "bad" dopplergrams with a manifestation of the linear interpolation in the primary 
data, the whole series was not processed. We have tested that the interpolation in the 
24-hours series strongly influences the behaviour of the k-u filter and the LOT program. 

In summary, 2.82 TB of primary data were processed. The result of the whole processing 
consists of 1004 horizontal full-disc two-component-velocity maps in Sanson-Flamsteed 
projection with an effective spatial resolution of 60". Each map has a size of 1520 x 1520 px, 
which makes together 18.5 GB of disk space. 
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Figure 19: Influence of the linear interpolation of missing frames in the series undergoing the 
temporal average oscillations removal. Left - dopplergram not influenced by the linear interpola- 
tion, right - p-modes of solar oscillation are clearly visible in the averaged dopplergram, which was 
calculated from the data series containing a significant amount of linearly interpolated frames. 
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Figure 20: The part of the power spectra of the dopplergram central area. The excess of power in 
the band of 0.32-0.34 px^^ in "bad" dopplergrams (dashed line) with respect to the "good" one 
(solid line) is clearly visible. 



5.2 Perspective removal 

The Sun is not a bright point and we observe it from the finite distance, therefore the 



positions are influenced by the perspective effect (Fig. 21). The point, which should be 
seen in the distance p from the axis, is due to the perspective in projection to the disc 
seen in the distance p + dr. Therefore the positions in the projection to the disc should 
be corrected by dr to get the correct positions in the heliographic coordinate system. The 



correction may be easily calculated from the geometry displayed in Fig. 21 
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From date - To date 


No. of day in series 


No. of processed series 


Z6 May lyyo — z4 Jui lyyo 


Dz 


1 1 Q 

iio 


io Apr iyy/ - 14 Jul iyy/ 


no 

yz 


1 O 1 

ioi 


y Jan iyyo — iU Apr iyyo 


yi 


1 7/1 


io Mar iyyy — Zo May iyyy 


IV 


yo 


97 Ma\r 9nnn — ThI 9nnn 




uo 


28 Feb 2001 - 28 May 2001 


88 


105 


18 Mar 2002 - 3 Jun 2002 


77 


83 


8 Aug 2003 - 23 Sep 2003 


46 


2 


18 Oct 2003 - 24 Nov 2003 


37 


22 


4 Jul 2004 - 6 Sep 2004 


64 


48 


25 Jun 2005 - 31 Aug 2005 


67 


86 


24 Mar 2006 - 22 May 2006 


59 


22 


Summary 


806 


1004 



Table 1: The overview of the processed data. Note that on each suitable day, two 24-hours data 
series sampled by 12 hours were processed. 



Two equations have to be solved: 

dr + p dr 

— ^=tan(ZJ=— , (25) 

d w 

and 

R'^ = p'^ + w^. (26) 

As the result we get 

We assume that w is always positive, because only one hemisphere of the Sun is relevant 
for the observations from the Earth or Li point, where SoHO is located. 

All the images have to be corrected to this effect, otherwise the displacement due to the 
perspective effect may reach the value of 1 px, which corresponds to 0.12° in heliographic 
coordinates. 



5.3 Oscillations removal 

From all frames, the line-of-sight component of the Carrington rotation is subtracted and 
the effect of a perspective is corrected. The frames are transformed so that the heliographic 
latitude of the disc centre 6o = and the position angle of the solar rotation axis P = 0. 
In the dopplergrams, the signal of p-modes of the solar oscillations is very well visible. The 
pattern of oscillations causes a lot of troubles during the application of the local correlation 
tracking method. To exclude this, the p-modes are suppressed using weighted average (see 
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Figure 21: Perspective effect. 



Hathaway, 1988). The weights have a Gaussian form given by the formula (16). We sample 
averaged images in the interval of 15 minutes. The filter suppresses more than five hundred 
times the solar oscillations in the 2-4 mHz frequency band. 

The frames with suppressed signal of solar oscillations are processed in the way de- 
scribed in detail in Section |3j Here I bring just a brief summary. 

The processing of averaged frames consists of two main steps. In the first main step the 
mean zonal velocities are calculated and, on the basis of expansion to the Fay's formula 
u = cq + ci sin^ b+C2 sin^ b, the differential rotation is removed. In the second main step, the 
LCT algorithm with an enhanced sensitivity is applied. Finally, the differential rotation 
(obtained in the first step) is added to the vector velocity field obtained in the second main 
step. Both main steps can be divided into several sub-steps, which are mostly common. 

1. The data series containing 96 averaged frames is "derotated" using the Carrington 
rotation rate in the first step and using the calculated differential rotation in the 
second step. 

2. Derotated data are transformed into the Sanson-Flamsteed coordinate system to 
remove the geometrical distortion caused by the projection to the disc. The Sanson- 
Flamsteed (also known as sinusoidal) pseudo-cylindrical projection conserves the ar- 
eas and therefore is suitable for the preparation of the data used by LCT. 



3. Remapped data undergo the k-u filtering (e. g. Title et al. , 1989) with the cut-off 
velocity 1 500 m s~^ for suppression of the noise coming from the evolutionary changes 
of supergranules, of the numerical noise, and for the partial removal of the "blind 
spot" (an effect at the centre of the disc, where the supergranular structures are 
almost invisible in dopplergrams due to the prevailing horizontality of their internal 
velocity field). 
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Finally, the LCT is applied: the lag between correlated frames is 4 hours, the cor- 
relation window has FWHM 60", the measure of correlation is the sum of absolute 
differences and the nine-point method for calculation of the subpixel value of dis- 
placement is used. The calculated velocity field is averaged over the period of one 
day. 



5. The resulting velocity field is corrected using the formula (15) in Section 3.3 The 
directions of the vectors before and after the correction are without change. Finally, 
fj, component is corrected for the data-processing bias of —15 ms~^ determined in 
Section 13.31 



5.4 Visualization 



The processed data may be visualised in many ways. Some of them that were used during 
the analysis of the data are shown in this subsection. As an example the 24-hours averaged 
horizontal large-scale flow field centered at 00.15 UT of March 6th 2001 is taken. 



The commonly used type of visualization is the vector arrow-plot (Fig. 22). This type 



of visualization is not useful for the study of the details at the full disc, but provides an 
important overview of the behaviour of the detected flows. The same field with mean flows 
(differential rotation, meridional circulation) removed is in local helioseismology usually 
called as the "sub-surface weather" (SSW). Since motions of the supergranules are mea- 
sured on the surface of the Sun, it is probably possible to call it a "surface weather". It 
provides an overview about the motions with respect to the very-large-scale background. 
An additional information comes from the map of absolute magnitudes of the velocities 



(Fig. 23). 



-0. 






.0 

-50.0 ' 
-100,0 
-1 50.0 
-200,0 
-250.0 
-300,0 
-350.0 
-400,0 



Figure 
March 
posed. 



22: Left - Arrow-plot of the surface horizontal large-scale flow field calculated for 6th 
2001. Right - '^Surface weather^^ at the same instant. The MDI magnetogram is superim- 
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Figure 23: Magnitude of the flow field. Note the increased speed in the northern (top) hemisphere 
near the central meridian. This area corresponds to the preceding polarity of the complex sunspot 
group. It shows that the preceding polarity of young sunspot group is moving faster than the 
following one and stretches the bipolar magnetic field. 

The largest flows consist of differential rotation and meridional circulation. The differ- 
ential rotation is calculated as the integral of the zonal component over the longitudes. To 
get a smooth curve, a parabolic fit in even powers of sinus of latitude - Eq. (|3]) - is usually 
calculated. This fit is performed mostly from the historical reason. As it was already 
mentioned in Section [2| this type of fitting has a lot of disadvantages. For example it does 
not introduce the asymmetries that can be present between the two hemispheres. 

The meridional circulation is calculated as an integral of the meridional component of 
the velocity field over longitudes. It brings some information about the transfer of the 
angular momentum and the magnetic flux towards the solar poles. It is considered mostly 
poleward with one cell per hemisphere. 



The mean flows calculated from the example result are displayed in Fig. 24 



The LCT technique basically provides only two components of the generally three- 
dimensional vector field. Nonetheless, it is possible to recover some information about the 
third (vertical) components also from the two-dimensional data. For this, the horizontal 



divergence Dh (Fig. 25 left) is calculated. 



Dh = - 



'dv^ ^ dvy 
dx dy 



Assuming the continuity equation in the form of V • v = we obtain: 

dvz 
dz 



{21 



(29) 
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Differential rotation fit: A=13.21. B = -0.23, C = -0.31 



100 [ 



Meridional flow fit: 5hift= 1 .82, slope= 1 .03 
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Figure 24-' Mean componets of the surface flow field. Left - differential rotation and its parabolic 
fit, right ~ meridional circulation and its linear fit. 



Therefore in the regions, where Dh is negative, we may expect a sink - i. e. a downflow 
and vice versa. 



An information about the twist comes from vorticity (Fig. 25 right). Vorticity Vz is 
defined as 

\4 = (V X v), . (30) 

It denotes the twist in the flow field. If Vz is negative, then the clockwise motion may be 
expected. Vortical structures are often present in areas of upflows and downflows. They 
are a consequence of the influence of the Coriollis force. 

The LCT technique was applied to the full-disc dopplergrams. From the resulting 
horizontal flow field (and neglecting the vertical component) it is possible to calculate the 
artificial dopplegram. 

Vd = Vh sin p, (31) 

where Vd is the line-of-sight velocity, Vh is the horizontal velocity, and p the heliocentric 
angle. 



Comparing the artificial dopplergram to the measured one (Fig. 26 ) provides an infor- 



mation about the difference of the plasma speed (measured by the line-of-sight velocity 
in measured dopplergrams) with the speed derived from the motion of the super granules. 



The difference (white and black areas near the disc edges in Fig. 26) may be caused by 
the wave-like motions or by the bias in the supergranules tracking procedure. This effect 
is known as the super-rotation of supergranules and was mentioned by several papers, first 



by Duvall (1980), lately confirmed by e. g. Snodgrass & Ulrich (1990), or recently by Beck 



& Schou (2000). The measured rotation of supergranules seems to be higher than the 
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Figure 25: Left - horizontal divergence of the flow field, right - vorticity of the flow field. 



rotation speed of the plasma in them. Although this super-rotation was recently explained 



as the effect of projection in the dopplergrams (Hathaway et al. , 2006), some uncertainties 




, 100,0 



r 




Figure 26: Maps of line-of-sight components. Left - artificial dopplergram constructed from the 
horizontal flow field with rigid rotation removed, right - original dopplergram with rigid rotation 
removed. 
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about this phenomenon still remain (e. g. Meunier & Roudier, 2007). When analysing the 



velocity powerspectra, some evidence about the traveling waves was found by |Gizon et al. 
(|2003l). 
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6 Long-term behaviour 

The subject of this chapter is a verification of the performance of the method described 
in Section |3] on the real data and the investigation of long-term properties of the fiows at 
largest scales obtained with this method. We shall also discuss the infiuence of magnetic 
fields on the measured zonal fiow in the equatorial region. 

The data are processed as described in Section [5j In this chapter, the interest was 
focused on the properties of the mean zonal and meridional components. Therefore, from 
each two-component horizontal velocity field the mean zonal and meridional components 
depending only on heliographic latitude were calculated as the longitudinal average of 
the flow map, using 135 longitudinal degrees around the central meridian. As stated in 
Section [s} the accuracy for each velocity vector is 15 ms~^ for velocities under 100 ms~^ 
and 25 ms~^ for velocities above 100 ms~^. These inaccuracies have a character of a 
random error, therefore for the mean zonal and meridional components the accuracy is in 
the worst case 1 ms""*^. 

The calculated surface flows may be biased by the projection effects, although the tests 



on the synthetic data did not show any signs of them. Hathaway et al. (2006) showed 
that the apparent superrotation of the structures tracked in dopplergrams reported by 
many studies can be explained as the projection effect. However, this bias would produce 
the systematic error, which should influence neither the periodic analysis, nor the relative 
motions of the active regions with respect to their surroundings. 



6.1 Results 

6.1.1 Long-term properties 

For the study of long-term evolution of surface flows maps, the mean zonal and merid- 
ional components were calculated. The maps of mean meridional component in time and 
heliographic latitude are shown in Fig. |27j It can be clearly seen that on the northern 
hemisphere the flow towards the northern pole dominates while on the southern hemi- 
sphere the flow towards southern pole prevails. The "zero line" , the boundary between the 
flow polarities, is not located exactly on the solar equator and seems to be shifted to the 



south in the period of increased solar activity (2001 and 2002). In agreement with Meunier 



et al. (1997) and Cameron & Hopkins (1998) the meridional flow is found stronger in the 



periods of increased solar activity by about 10 ms~^ than in periods with lower magnetic 
activity. 

A similar map was made for the zonal component. The mean equatorial zonal velocity 
for all the data is 1900 ms~^. For all the processed data the dependence on latitude is 
close to a parabolic shape, parameters of which change slowly in time. The residua of the 

★ This chapter was submitted to Astronomy & Astrophysics as Svanda, M., Klvaha, M., Sobotka, M., and 
Bumba, V., Large-scale horizontal flows in the solar photosphere. II. Long-term behaviour and magnetic 
activity response. 
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Figure 27: Mean meridional flow in time and heliographic latitude. It can be clearly seen that 
for almost all the processed measurements a simple model of one meridional cell per hemisphere 
would be sufficient. However, some local corruptions of this simple idea can be noticed on both 
hemispheres. 
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Figure 28: Torsional oscillations. The residua of the mean zonal flow with respect to its parabolic 
fit displayed in time and heliographic latitude. In the period of weak magnetic activity the pattern 
of belts propagating towards the equator is very clear. In the periods of stronger magnetic activity 
the flow field is influenced by local motions in active regions and therefore the pattern of torsional 
oscillations is not clearly seen. 
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zonal velocity with respect to its parabolic fit given by 

ff, = + ai6 + a2&^, (32) 

where b is the heliographic latitude and vj, the mean zonal velocity in the given latitude, 
were calculated in order to see if the torsional oscillations can be detected in the measure- 



ments. As it is displayed in Fig. |28} the method clearly reveals torsional oscillations as an 
excess of the mean zonal velocity with respect to the zonal velocity in the neighbourhood. 
The behaviour of torsional oscillations is in agreement with their usual description - the 
excess in magnitude is of the order of 10 ms~^, they start at the beginning of the solar cycle 
in high latitudes and propagate towards the equator with the progression of the 11-year 
cycle. However, due to the used method, the visibility of torsional oscillations decrease 
with increasing solar activity. In the periods of strong activity both belts are not so clearly 
visible since the large-scale velocity field and its parabolic fit are strongly influenced by the 
presence of magnetic regions. However, the torsional oscillations belts still remain visible 
when the mean zonal component is symmetrised with respect to the solar equator. The 
meridional flow or torsional oscillations depending on time and latitude are not in the focus 
of this study, they are just used to check the ability and performance of our method. 

6.1.2 Periods in the mean components 

The mean zonal and meridional components in the equatorial area (averaged in the belt 
b = —5° — 1-5°) were analysed in order to examine the periods contained in the data. 
Since the data are not equidistant at all, a simple harmonic analysis cannot be used, so 



that the Stellingwerf method (Stellingwerf , 1978) was applied. It works on the principle of 



the phase dispersion minimization. The method sorts the data for every searched period 
into the phase diagram. Then the phase diagram is divided in a few (mostly ten) parts 
and for every part the mean dispersion is calculated. If the studied period is reasonable, 
the data-points group along the periodic curve and the dispersion in each part of the phase 
diagram is smaller than the dispersion of the whole data series. The normalized parameter 
9 G (0, 1) describes the quahty of a given period. 

The influence of the calculated flow fields by the the position and orientation of the solar 
disc (position angle of the rotation axis P, heliographic latitude of the centre bo) cannot 
be completely excluded, so the series of both parameters sampled on the same dates when 
the measurements of surface flows exist also passed the period analysis. Also the periods 
caused by sampling of the data using the same method were verified. Periodograms are 



displayed in Fig. 29 



Unfortunately, it has to be concluded that no significant period in the available data 



set was detected. As it can be seen in Fig. 29, there exist non- convincing (the values of 



the parameter 6 are quite high, which means that probably the periods are not significant) 
signs of periods detected in the real data, which are not present in the control data set (P, 
bo). Values of the suspicious periods are 657 days (1.80 years) in the meridional component 
and 1712 days (4.69 years) in the zonal component. Note that the 1.8-year period was also 



detected by Knaack et al. (2004). It is claimed to be related to a possible Rossby wave 
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Figure 29: Periodograms determined using the Stellingwerf method. Parameter 9 signifies the 
normalized phase variation. Upper left: Periodogram of mean equatorial zonal velocity. Upper 
right: Periodogram of mean equatorial meridional velocity. Bottom left: Periodogram of sampled 
heliographic latitude of the centre of the solar disc. Bottom right: Periodogram related to the 
sampling of data. 



r-mode signature in the photosphere with azimuthal order m ~ 50 reported by |Kuhn et ah 
(2000), but lately disputed e.g. by Williams et al. (2006). The period estimate for such 
an r-mode is close to 1.8 years. According to Knaack et al. (2005), such a periodicity 



was observed in the total magnetic flux only on the southern hemisphere from 1997 to 
2003. The coupling between the zonal flow and the meridional circulation could transfer 
the signal of the r-mode motion to the mean meridional component. 

The detected suspicious periods may not be of solar origin. The sparse data set suffers 
from aliasing caused by a bad coverage of the studied interval. To confirm the periods, a 
far more homogeneous data set is needed. This may be a task for the ongoing space borne 
experiment Helioseismic Michelson Imager (HMI), which will be a successor of MDI. 

Detected periodicities are absent in Mt. Wilson torsional oscillations time series, which 



is a far more homogeneous material than the one used in this study. Snodgrass fc Ulrich 



(1990) also did not find any time variations in the study tracking the features in the low 



resolution dopplergrams covering homogeneously 20 years of Mt. Wilson observations. All 
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Figure 30: Mean zonal equatorial velocity versus the sunspot area in the near- equatorial belt. 
We decided to divide the data in two regimes along the velocity axis. Although the division is 
arbitrary, we believe that it is supported by the theory of the dynamical disconnection of sunspots 
from their roots. 

the arguments written above led to leave the detected periods as suspicious, as they cannot 
be confirmed from the current data set. 



6.2 Relation to the magnetic activity 

The coupling of equatorial zonal velocity (average equatorial solar rotation) and the solar 
activity in the near-equatorial area (belt of heliographic latitudes from — 10 ° to -1-10 °) was 
also investigated. The average equatorial zonal velocity incorporates the average super- 
granular network rotation and also the movement of degenerated supergranules influenced 
by a local magnetic field with respect to their non-magnetic vicinity. Indexes of the solar 
activity were extracted from the daily reports made by Space Environment Center National 
Oceanic and Atmospheric Administration (SEC NCAA). Only the days when the measure- 
ments of horizontal flows exist were taken into account. As the index of the activity we 
have considered the total area of sunspots in the near-equatorial belt and also their type. 

First of all, the correlation coefficient p between the mean equatorial zonal velocity 
and the sunspot area in the near-equatorial belt was computed with a value of p = —0.17. 
There is no significant linear relation between these two indices. The dependence of both 



quantities is plotted in Fig. 30 There can be clearly found two different regimes, which 
are divided by the velocity of approximately 1890 ms~^. In one regime (77 % of the cases), 
the equatorial belt rotates about 60 ms~^ faster (1910 ±9 ms~^; hereafter a "fast group") 
than the Carrington rotation, in the other one (23 %) the rotation rate is scattered around 
the Carrington rate (1860 ±20 ms~^; hereafter a "scattered group"). The division in these 
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Figure 31: Sunspot area sampled in the same times when the measurements of the horizontal 
flows exist. Two regimes of the near- equatorial belt rotation are displayed. Diamonds denote the 
"fast" rotating equatorial belts, crosses the "scattered" group. 



two suggested groups using the speed criterion is arbitrary. If there exist only two groups, 
they certainly overlap and only a very detailed study could resolve their members. One 
may also see more than two groups in Fig. 30, The arguments for the division in just two 
groups will follow. 

For both regimes a typical sunspot area does not exist. The distribution of both regimes 
in time is displayed in Fig. 31 The data in the periods of larger solar activity (years 2001 
and 2002, these are also the only years when the data cover two Carrington rotations 
continuously) show that both regimes alternate with a period of one Carrington rotation. 

The histogram of the mean zonal equatorial velocity has a similar, i. e. bimodal, char- 
acter like in Fig. 30 with a greater second peak, because such a histogram is constructed 
not only from belts containing magnetic activity but also from the belts where no mag- 
netic activity was detected. The mean equatorial rotation for all the data is 1900 ms~^, 
1896 ms~^ for the equatorial rotation in the presence of sunspots, and 1904 ms~^ for days 
without sunspots in the equatorial region. 

Such bimodal velocity distribution is in disagreement with the results obtained by time- 
distance helioseismology by Zhao et al. (2004). In this work the authors found that the 
stronger the magnetic field the faster such a magnetic element rotates. They observed about 
70 ms~^ faster rotation than the average for magnetic areas with magnetic fields stronger 
than 600 G. It might be possible that the size of the magnetic area and its magnetic field 



strength play different roles in the influencing of the plasma motions. However, Meunier 



(2005) showed that the size of the magnetic area and the maximum magnetic field strength 
or the total magnetic flux correlate quite well. 

Detailed studies of the sunspot drawings obtained from the Patrol Service of Ondfejov 
Observatory and the Mt. Wilson Observatory drawings archive revealed that in the "fast" 
group, the new or growing young active regions were present in the equatorial belt. On 
the contrary, in the "scattered" group the decaying or recurrent active regions prevailed in 
the equatorial area. The deceleration of the sunspot group with its evolution was noticed 
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Figure 32: Left: Distribution of two equatorial rotation modes in the year 2002. Right: Deriva- 
tives of the mean zonal velocity (solid curve) and the sunspot area in the near-equatorial region 
(dashed curve) in 2002. Both quantities correlate with each other quite nicely. 




Figure 33: Active region morphological types ('left ) and the number of active regions in the near- 
equatorial belt distributed as a function of time and sunspot area ( rightj. 



e.g. by Ruzdjak et al. (2004). Moreover, the results suggest that the new and rapidly 
growing sunspots in the studied sample (March to May 2001 and April to June 2002) 
move with the same velocity. This behaviour could be explained by an emergence of 
the local magnetic field from a confined subphotospheric layer. According to the rough 
estimate (Corbard & Thompson, 2002) the speed of 1910 it 9 ms~^ corresponds to the 
layer of 0.946 ±0.008 Rq, where the angular velocity of rotation suddenly changes. During 
the evolution, the magnetic field is disrupted by the convective motions. An interesting 
behaviour is displayed by the alternation of the "fast" and "scattered" regimes (see Fig. 32 ) 
with the period of one Carrington rotation. It suggests that active regions in the equatorial 
region emerge in groups. We have to keep in mind that this study produces a very rough 
information due to the averaging of all effects in the equatorial belt. 

The observed behaviour could be a manifestation of the disconnection of magnetic field 
lines from the base of the surface shear during the evolution of the growing sunspot group. 
This behaviour was theoretically studied by Schiissler & Rempel (2005). They suggested 
the dynamical disconnection of bipolar sunspot groups from their magnetic roots deep in 
the convection zone by upflow motions within three days after the emergence of the new 
sunspot group. The motion of sunspots changes during those three days from "active" to 
"passive" . The active mode is displayed by motions reasonable faster with respect to the 
non-magnetic surroundings. The passive mode means mostly the deceleration of sunspot 
motions and influence of the sunspot motions only by the shallow surface plasma dynamics. 
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The theory of the disconnection of sunspot groups from their magnetic root supports the 
division of the data set in two groups. 



As an example, the active region NOAA 9368 (Fig. 34) was selected to show the be- 
haviour of large-scale velocities in time. One see that the leading part of the active region 
rotates faster than the surroundings in the first day of observation and the whole group 
slows down in next two days. The inspection of details in the behaviour of selected active 
regions in the whole data set will be the subject of the ongoing studies. 

The division of the equatorial belt into 10 sectors helped the investigation of the be- 
haviour of different active region types. The dependence of the mean zonal velocity in 
the sector containing the studied active region on the morphological type of active region 
is shown in Table |2} It can be clearly seen that, in average, more evolved active regions 
rotate slower than less evolved or young active regions, what is in agreement with e. g. 



Balthasar et al. (1986). No particular behaviour for various active region types was found 



(see Fig. 33). 

We have also studied how the presence of the magnetic active areas will influence the 
average flow field. Since we found that a direct correlation is weak due to the existence of 
two different regimes, we decided to study the temporal change of both quantities. The aim 
is to study whether an emerging active region in the near-equatorial belt will influence the 
average equatorial rotation. We computed numerical derivatives of the total sunspot area 
in the near-equatorial belt and of the average zonal equatorial flow. We have found that 
the correlation coefficient between both data series is p = 0.36 and is higher for the "fast 
group" (p = 0.41) than for the "scattered group" (p = 0.24). The correlation is higher in 
periods of increased magnetic activity in the equatorial belt. For example, for data in the 
year 2001 the correlation coefficient is P2001 = 0.58 and for the year 2002 P2002 = 0.52; see 



Fig. 32, In both particular cases, the correlation is higher for the "fast regime" (p ~ 0.7) 



than for the second group. 



Table 2: Mean synodic rotation velocities of different active regions types and the average equato- 
rial rotation of all data. The measured speed values may be systematically biased by the projection 



effect (see Hathaway et al., 2006). 



Sunspot 
type 



Mean rotation 

[ms^-*-! 



A 


1^ 


^93 ± 48 


B 


1^ 


^93 ± 49 


C 


1^ 


^90 ± 71 


D 


1^ 


^80 ± 73 


E 


1^ 


^80 ± 50 


F 


1^ 


^74 ± 73 


H 


1^ 


^72 ± 51 


Average 




1900 
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Figure 34-' Case study: Evolution of the flows in and around the active region NOAA 9368 on 
March 6, March 1, and March 8, 2001. The leading polarity rotate significantly faster than the 
following one and the non-magnetic surroundings in the first day. The whole group slows down in 
the other two days. As the background image, the MDI magnetogram smoothed to the resolution 
of the measured flow field is used. The magnetic field intensities are displayed in the range of 
— 1800 (black) to +1800 (white) Gauss in the linear scale. 
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It is important to know that LCT in our method measures basically the motions of 
supergranules influenced by the magnetic field, not by spots that are recorded in the used 
solar activity index. Therefore the results can be biased by the fact that the presence of 
magnetic field does not necessarily mean the presence of a sunspot in the photosphere. 

We think that despite an apparent disagreement, our apparently conflicting results can 
be valid and are in agreement with the results published earlier. Basically, as described 
by e.g. Hathaway & Wilson (1990) and explained in the model of Brun (2004), the 



solar rotation in the lower latitudes is slower in the presence of magnetic field 
summarized in Table |2] and displayed in Fig 



30 



This is 

In most cases, "spotty" equatorial belts 
seem to rotate slower than the average for the whole data series. However, it is clear that 
emerging active regions cause in most cases the increase of the rotation rate. This is in 



agreement with a generally accepted statement found first by Howard & Harvey (1970) 



and Golub & Vaiana (1978). The relation, obtained using a linear fit on our data set, can 



be described by the equation 



Av ~ 0.2 AAsunspots ms" 



(33) 



where Av is a change of the equatorial rotation speed with respect to the Carrington 
rotation and AAsunspots is a change of sunspot area in the equatorial belt (in 10~^ of solar 
hemisphere). We estimate that strong local magnetic areas rotate few tens of ms~^ faster 
than the non-magnetic surroundings. 

The results seem to be in contradiction to the previous study by [Snodgrass &: Ul- 



rich (1990), who tracked the features in the low- resolution dopplergrams measured at the 



Mt. Wilson 46 m tower telescope for the period of 20 years. The authors interpret the 
detected flows as the velocity of the supergranular network, although the spatial resolution 
was lower than the size of individual supergranules. They found two regimes taking place 
in the rotation of the photosphere - the quiet Sun and the active regions. The results 
showed that the regions occupied by the magnetic field displayed a slower rotation than 
the non-magnetic vicinity. Title et al. ( [I992j ) used the high resolution data obtained at 
the Swedish Vacuum Solar Telescope on La Palma, Canary Islands, and found that the 
velocity field on the granular scale is different in the regions of the quiet Sun and in an 
active-region plage. Horizontal flow speeds measured by LCT were a factor of 2 slower in 
magnetic field regions than in the quiet Sun. 



6.3 Conclusions 

We have verified that the method developed and tested using the synthetic data is suitable 
for application to real data obtained by the MDI onboard SoHO and maybe also to the 
data that will be produced by its successor Helioseismic Michelson Imager (HMI) onboard 
the Solar Dynamic Observatory (SDO). HMI will have a greater resolution and will cover 
larger time span than two months each year. We verified that the long-term evolution of 
the horizontal velocity fields measured using our method is in agreement with generally 
accepted properties. 
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During the periodic analysis of the equatorial area we found two suspicious periods in 
the real data, which are not present in the control data set containing the inclination of the 
solar axis towards the observer, the quantity that can bias systematically and periodically 
the results by a few ms~^. The periods of 1.8 year and 4.7 years need to be confirmed 
using a more homogeneous data set. 

We also found that the presence of the local magnetic field generally speeds-up the 
region occupied by the magnetic field. However, we cannot conclude that there exists a 
dependence of this behaviour for different types of sunspots. We can generally say that 
the more evolved types of active regions rotate slower than the young ones, however the 
variance of the typical rotation rate is much larger than the differences between the rates 
for each type. We have found that the distribution of active regions rotation is bimodal. 
The faster-rotating cases correspond to new and growing active regions. Their almost 
constant rotation speed suggests that they emerge from the base of the surface radial shear 
at 0.95 Rq. The decaying and recurrent regions rotate slower with a wider scatter in their 
velocities. This behaviour suggests that during the sunspot evolution, sunspots loose the 
connection to their magnetic roots. Both regimes alternate with a period of approximately 
one Carrington rotation in years 2001 and 2002, which suggests that new active regions 
emerge in groups and may have a linked evolution. 
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7 MOTIONS AROUND THE FOOTPOINTS OF THE ERUPTIVE FILAMENT 



7 Motions around the footpoints of the eruptive filament 



Dynamic processes on the Sun are linked to the evolution of the magnetic field as it is 
passing through the different layers from the convection zone to the solar atmosphere. 
In the photosphere, magnetic fields are subject to diffusion due to supergranular flows 
and to the large-scale motions of differential rotation and meridional circulation. The 
action of these surface motions on magnetic flelds plays an important role in the formation 
of large-scale fllaments (Mackay & Gaizauskas, 2003). In particular, the magnetic flelds 



that are transported across the solar surface can be sheared by dynamic surface motions, 
which result in shearing of the coronal fleld. This corresponds to the formation of coronal 



flux ropes in models which can be compared with Ha fllament observations (Mackay & 



van Ballegooijen, 2006b). Many theoretical models try to reproduce the basic structure 



and the stability of fllaments by taking into account surface motions as quoted above. 
These models predict that magnetic flux ropes involved in solar fllament formation may be 
stable for many days and then suddenly become unstable resulting in a fllament eruption. 
Observations show that twisting motions are a very common characteristic of eruptive 
prominences (see for example Patsourakos &; Vial, 2002). However, it is still unknown 



whether the magnetic flux ropes emerge already twisted or if it is only the photospheric 
motion that drives the twisting of the fllament magnetic fleld. Therefore, the mechanisms 
that drive fllament disappearance remain uncertain. The destabilization can come from 
the interior of the structure or by an outside flare. 



In the paper by Rondi et al. (2007), local horizontal photospheric flows were measured 



at high spatial resolution (0.5") in the vicinity of and beneath a fllament before and during 
the fllament's eruptive phases (the international JOP178 campaign). It was shown that 
the disappearance of the fllament initiates in a fllament gap. Both parasitic and normal 
magnetic polarities were continuously swept into the gap by the diverging supergranular 
flow. We also observed the interaction of opposite polarities in the same region, which could 
be a candidate for initiating the destabilization of the fllament by causing a reorganization 
of the magnetic fleld. Here we investigate the large-scale photospheric flows at moderate 
spatial resolution (2") beneath and in the vicinity of the same eruptive fllament. 



7.1 Observations 

During three consecutive days of the JOP 178 campaign, Oct 6, 7, and 8, 2004 (http:/ /gaia. 
bagn.obs-mip.fr/jopl78/index.html), we observed the evolution of a fllament that was close 
to the central meridian. We also observed the photospheric flows directly below the fllament 
and in its immediate area. The filament extends from —5° to —30° in latitude. A filament 
eruption was observed on October 7, 2004 at 16:30 UT at multiple wavelengths from 

★ This chapter was done together with Thierry Roudier from Laboratoire d'Astrophysique de I'Observatoire 
Midi-Pyrenees, Tarbes, France, and is a part of a paper submitted to Astronomy & Astrophysics as 
Roudier., Th., Svanda, M., Meunier, N., Keil, S., Rieutord, M., Malherbe, J. M., Rondi, S., Molodij, G., 
Bommier, V., and Schmieder, B., Large-scale horizontal flows in the solar photosphere. III. Effects on 
filament destabilization. 
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MDI 9h40 UT 




Longitude 

Figure 35: MDI longitudinal magnetic field data on October 7, 2004 beginning of the 

sequence at 9:40 UT and at its end at 22:40 UT. In order to see the magnetic field evolution we 
overploted the shape of the filament at 13:30 UT in both figures. 



ground and space instruments. The eruption produced a Coronal Mass Ejection (CME) at 
approximately 19:00 UT that was observed with LASC0-2/S0H0 and two ribbon flares 
observed with SOHO/EIT. MDI/SOHO longitudinal magnetic field and Doppler velocity 
were recorded with a cadence of one minute during the 3 days. The Air Force ISOON 
telescope located at the National Solar Observatory/Sacramento Peak provided a full-disc 
Ha image every minute. The pixel sizes were respectively 1.96" for MDI magnetograms 
and dopplergrams and 1.077" for ISOON Ha images. 

Our primary goal was to derive the horizontal fiow field below and around the filament. 
Co-alignment between SOHO/MDI magnetograms and ISOON data was accomplished 
by adjusting the chromospheric network visible in Ha (ISOON) and the amplitude of 
longitudinal MDI magnetograms to an accuracy of one pixel (1.96"). The general magnetic 



context before and after the filament eruption is shown Fig. 35 



The local correlation tracking technique was applied to the set of full-disc MDI dopp- 
lergrams measured in three days around the time of filament eruption in high cadence. So 
the data in this campaign are suitable to be processed with the procedure described in 
Section |5| Due to the time series length and the required temporal resolution we had to 
adjust some parameters of data processing. The most significant difference is a sampling 
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Figure 36: Horizontal photospheric flow field derived using LCT applied to full-disc dopplergrams. 
The filament observed by ISOON on October 7, 2004 ^-^ superimposed. The hexagon indicates the 
location where the filament eruption started. 



of averaged p- modes-free frames in 1-minute cadence instead of 15 minutes. The averaged 
frames were tracked using the Carrington rotation rate (with an angular velocity of 13.2 
degrees per day), so that all the frames have the same heliographic longitude of the central 
meridian [Iq = 62.24°). The LCT method applied to full-disc dopplergrams is character- 
ized by a Gaussian correlation window {FHWM = 60") and a time lag between correlated 
frames of 1 hour (basically 60 frames). In all cases, one half of the intervals were before 
the eruption and the second half after the eruption. All the pairs of correlated frames in 
the studied intervals were averaged to increase the signal-to-numerical-noise ratio. 



7.2 Photospheric flow pattern below and around the filament 

In this subsection we describe the flows associated with the filament eruption, with particu- 
lar emphasis on the filament evolution and the mean east-west (zonal) velocities properties. 



The 13- hours averaged flow fields in the vicinity of the filament are shown in Fig. 36 
The measurements were compared with similar measurements based on tracking of indi- 
vidual magnetic features in MDI magnetograms. All the methods provide similar results. 
The large-scale flows are well structured and show both converging and diverging velocity 
patterns. We observe in particular a large-scale stream in the north-south direction par- 
allel to the filament located about 10° to the east between —20° and —30° in latitude and 



58° and 47° in longitude. This flow stream is clearly visible in Fig. 36, and its dynamics 
can be seen at http:/ /gaia. bagn.obs-mip.fr/jopl78/oct7/mdi/7oct-mdi. htm. Around —20° 
in latitude the velocities of the differential rotation amplitude start to dominate. How- 
ever, in both measurements we observe that the north-south large-scale stream on the 
eastern edge of the filament disturbs the regular differential rotation. The location of the 
north-south stream is close to the location where the filament eruption begins (longitude 
/ = 56°, latitude b = —26° in Carrington coordinates). The amplitude of the southward 
motions is approximately 40 ms~^, which is of the same order as the mean observed flows. 
Below the latitude of —20°, the combination of differential rotation and the north-south 
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Figure 37: Evolution of the filament during its eruption around 16:30 UT on October 7, 2004- 
The arrow indicates a fixed position for all the subframes. A and B denote two parts of the 
filament. 



stream cause opposite polarities to move closer, which strongly increases the tension in the 
magnetic field very close to the starting point of the filament eruption. 



The filament's evolution can be seen in Fig. 37 The north-south stream fiow visible in 



Fig. 39 (left) crosses over the part of the filament labeled A. The arrow on Fig. 37 indicates 



the same fixed point (325", 167") in all of the subframes. We observe a general southward 
motion of both the A and B segments of the filament. More precisely, we measure a tilt 
of these two filament segments at the point of their separation. Between 16:07 UT and 
16:58 UT the long axis of segment A rotates by an angle of 12° (clockwise) relatively to 
its western end, and the long axis of segment B of the filament rotates by an angle of 5.5° 
(clockwise) relatively to its western end. These rotations are compatible with the surface 



fiow shown in Fig. 39 (left) and in particular the north-south stream fiow. 

We constructed the mean zonal velocity profile, i. e. a profile of the differential rotation, 
to see if there is any correspondence to the presence of the filament. In this profile we 
found a strong secondary maximum visible at —23° of latitude. The secondary maximum 
indicates a decrease in the amplitude of the Vx component because the fiows in that region 
are oriented more in the north-south direction. That is partly due to the presence of the 
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Figure 38: Profiles for 7 Oct in different longitudinal belts: solid line and squares for —25° to 
— IT degrees with respect to the centre of the disc (not Carrington coordinates), dashed line and 
triangles for (F to 20°. 



north-south stream described above and to the local organization of the flow. In particular, 
converging and diverging flows in this region seem to have a north-south orientation. 

To distinguish the effects of the north-south stream from differential rotation, we com- 
puted the mean zonal velocities from —25° to —17° in longitude with respect to the cen- 
tre, where the north-south stream is present, and from 0° to 20° in longitude, where the 
differential rotation is not influenced by the current. The mean zonal velocities in the 
longitudinal belt, where the north-south stream is visible, exhibits clearly a secondary 
maximum (Fig. 38 ) indicating that the solar rotation rate at this location is closer to that 
of the equator. As a consequence, the plasma in the north-south stream, which transports 
magnetic structures, rotates faster at about —23° latitude, than do the magnetic structure 
located in the belt of longitude between 0° to 20°. The combination of these different sur- 
face motions (stream and differential rotation) tends to bring together fields with opposite 
polarities, this in turn constrains the magnetic field lines. 

We noted that the location of the starting-point of the filament eruption is around —26° 
in latitude, which is very close to that secondary maximum in the mean zonal velocity. 
Thus surface motions that bring together opposite polarities may play a role in triggering 
the filament eruption. 



7.3 Flow fields before and after the eruption 

In this subsection we discuss the properties of the flow field just before and after the 
filament eruption, which occurred at about 16:30 UT. At the point where the filament 
eruption begins (/ = 56°, b = —26° in Carrington coordinates), we detected a steepening 
of the gradient in the differential rotation curve. During the eruption, the gradient flattens 
out and a dip forms. Although the differential rotation curves describe the mean zonal 
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Figure 39: Horizontal motions measured before ('left ) and after ( right ) the eruption in the wide 
field of view. The filament observed by ISOON on October 7 2004 ^.^ 13:30 UT is superimposed. 
Bottom - The differences in the directions of both vector fields are evaluated using the weighted 



cosine defined by formula (24) in the sliding window with size of 3.5 heliographics degrees. It is 
seen that the worst correspondence of both vector fields is in the area of the north-south current. 
It demonstrates that the flow field in the field-of-view remained more-or-less stable during the 
filament eruption except for the close vicinity of the starting point. 
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Figure 40: The 3-hour average of the flow field in the close vicinity of the starting point (I = 56°, 
b = —26°) before (left) and after (^rightj filament eruption. The filament observed by ISOON on 
October 7 2004 15.20 UT is superimposed. The red boxes denote areas used for the zonal shear 
calculation. 



velocities on the full disc, the change of its gradient signifies the change in the stretch 
influencing the magnetic field in the loop over the area under study. We can express the 
surface rotation as an even power of sin 6 : lu = A + B sin^ b + C sin 6, where A is the 
angular velocity rate of the equatorial rotation and b is the heliographic latitude. From 
the data we find that the constant values (with their errors in parentheses) are before the 
eruption A = 13.375(0.010), B = -1.46(0.10), C = -1.42(0.20) and after the eruption 
A = 13.404(0.010), B = -1.78(0.10), C = -1.24(0.20). All the rates are synodical 
in degday"^. The full-disc profiles did not change significantly from before to after the 
eruption. For example, for a latitude of —30° the zonal velocity has values of 12.92 (resp. 
12.88) degday^^ (—34 ms~^, resp. —39 ms~^ in the Carrington coordinate system), for 
a latitude of —20° the values are 13.18 (—2 ms~^) resp. 13.18 (—3 ms~^) degday"^. 
Although the parameters of the smooth fitted curve did not change too much, the local 
residual with respect to the smooth curve changed at the latitude where the filament 
eruption starts. 



Fig. 39 displays the horizontal flows before and after the eruption in a wide field of view 
measured using the LCT method on supergranular structures and averaging the resulting 
velocities over 3 hours. Before the eruption, we can clearly see the north-south stream 
parallel to and about 10° to the east of the filament. This stream disturbs differential 
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Figure 41- Temporal evolution of the velocity shear in zonal components. The eruption of filament 
took place at 16:30 UT. 

rotation and brings plasma and magnetic structures to the south. Although differential 
rotation tends to spread the magnetic lines to the east, the observed north-south stream 
tends to shear the magnetic lines. After the eruption, only a northern segment of the 
filament is visible and the north-south stream has disappeared. 



Fig. 40 shows the flow field in more detail, at the site where the eruption starts. The 
shear in the zonal component at the point where the eruption starts (/ = 56°, h = —26° 
in Carrington coordinates) is clearly visible and prevails before and after the eruption, 
although the shape of the apparent vorticity has changed. This location corresponds to 
the area of upflows observed in the Meudon Ha dopplergram (Fig. 8 in Rondi et al. 



2007). 



We measured the shear as a difference between the mean zonal component in the 
areas just to the north and just to the south of the starting point. We measured the 
average zonal flow over boxes with a side of 2.3° located 2.9° to the north and to the 
south of the point where the filament eruption appeared to start. The evolution in the 
shear velocity measured as the difference between the mean flow in the two boxes as a 
function of time can be seen in Fig. 41 Six 2-hours averages of the flow fields were used to 



determine this figure. One can see that the shear velocity is increasing before the eruption 
and decreasing after the eruption. One hour before the eruption, the shear reached the 
value of (120±15) ms~^ in a distance of 5.2° (62 000 km in the photosphere). After the 
filament eruption, we observe a restoration of an ordinary differential rotation below —30° 
in latitude. 



7.4 Evolution on 6, 7 and 8 October 2004 



In Fig. 42 we compare the flow field in the filament area for the day of the eruption 
(October 7, 2004), with the flow fields on the preceding and following days. We see that 
the topology of the flows changed in the filament area: the daily evolution of the mean 
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Figure 4^-' Horizontal flow evolution in three days 6 ( upper leftj, 7 ( upper rightj and 8 ('bottom 
leftj October 2004- Bottom right - Mean zonal velocities, computed between —35^ and +20P in 
longitude and —3(F to (F in latitude. Dashed line (and triangles) represents mean zonal velocity 
for the LCT applied on dopplergrams on 6 Oct 2004- Solid line (and squares) represents mean 
zonal velocity for the LCT applied on Doppler on 7 Oct 2004- Dash-dot line (and diamonds) 
represents mean zonal velocity for the LCT applied on Doppler on 8 Oct 2004- 



zonal profiles can be seen in Fig. 42 (bottom right). The dashed line with triangles is 
the mean zonal velocity for 6 October 2004. This profile is relatively flat probably due to 
the short time sequence as only 3 hours of data were available. The differential rotation 
proflle for 7 October shows a secondary maximum around —23° in latitude, as discussed 
before. The 8 October proflle exhibits the same trend, but with a smaller amplitude and 
an eastward velocity for latitudes greater than —10°. The secondary maximum appears 
strongly reduced indicating a restoration to a more regular differential rotation pattern in 
that zone. 

One day before the eruption, the shear began to form at the site, where the filament 
eruption is triggered (/ = 56°, h = —26° in Carrington coordinates). The north-south 
stream is also visible. Both phenomena may store free energy in the coronal magnetic field 
configuration. The topology of the fiow and the stream is different a day after the filament 
eruption suggesting that after the disappearance of the southern part of the filament the 
conditions in the photosphere below the filament became more relaxed. This may suggest 
the mutual coupling of the photospheric fiow and the configuration of the coronal mag- 
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netic field. To confirm this idea, high-cadence high-resolution images and magnetograms 
covering the eruption time would be needed. 

7.5 Discussion and conclusion 

Filaments and prominences are important complex structures of the solar atmosphere 
because they are linked to CMEs, which can influence the Earth's atmosphere and near- 
space environment. Surface motions acting on pre-existing coronal fields play a critical role 
in the formation of filaments. They appear to reconfigure existing coronal fields, they can 
twist and stretch them and thereby deposit energy in the topology of the coronal magnetic 
field. Photospheric motions can also initiate coronal magnetic field disruption. Surface 



motions play an important role in formation Type B filaments (Tandberg-Hanssen, 1974), 
which are located between young and old dipoles and are long stable structure. This class 
of filaments requires surface motions to gradually reconfigure preexisting coronal fields 



e. g. DeVore & Antiochos, 2000) 



In a paper by Rondi et al. (2007) authors removed all the large-scale flows in order to 
focus on smaller-scale flow, such as mesogranulation and supergranulation. In this study 
we have retained the large-scale flows in order to study their influence on the triggering of 
a filament eruption. 

The eruption started around 16:30 UT at the latitude around —25° where the mea- 
surements of the horizontal flows based on dopplergram tracking show a modification of 
the slope in the differential rotation of the plasma. It seems to be a consequence of the 
presence of a north-south stream along the filament's position, which is easily measured 
by tracing structures in the dopplergrams. 

The observed north-south stream has an amplitude of 30-40 ms~^. In the sequence 
of Ha images that record the filament's evolution, the part of the filament, which is in 
the north-south stream, is rotated in a direction compatible with the flow direction of 
the stream. This behavior suggests that the foot-points of the filament are carried by the 
surface flows. The influence of the stream is strengthened by differential rotation. We 
should keep in mind that the filament extends from —5° to —30° in latitude and that 
the northern part of the filament is subjected to a faster rotation than the southern part. 
The stream, with a contribution of the differential rotation, causes the stretching of the 
coronal magnetic field in the filament and therefore contributes to destabilisation of the 
filament. The topology of the north-south stream changed after the filament eruption and 
the filament almost disappeared. 

We have measured the increase of the zonal shear at the site where the filaments 
eruption begins before the eruption. After the eruption, the shear suddenly decreased. 
This result suggests that the shear in the zonal component of the flow field is the most 
important component of the surface flow affecting the stability of the coronal magnetic 
field, and can lead to its eruption, which in turn can drive active phenomena such as 
ribbon flares and CMEs. This evolution of the shear in the flow field is probably related 
to the re-orientation by 70° (or 110°) of the transverse field after the eruption seen in the 



daily vector magnetograms obtained with THEMIS (Rondi et al. , 2007). 
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All of the features observed in the topology of the horizontal velocity fields at the 
starting-point site could contribute to destabilise the filament resulting in its eruption. We 
propose that the stability and evolution of filaments are influenced by surface flows that 
carry the footpoints of the filament. 

Filaments or prominences are important complex structures of the solar atmosphere. 
Several mechanisms are probably involved in the creation of filament eruptions: the action 



of surface motions to create or increase the helicity of the flux rope (van Driel-Gesztelyi 
2005 or Romano et al. , 2005), reconnecting field lines in the corona (Mackay & van Bal- 



legooijen, 2006a), the chirality evolution of the barb (Su et al. 2005), etc. The coronal 



magnetic field is generally thought to be anchored in the photosphere and the flux trans- 



port on the solar surface (Wang et al. , 1989) is the natural mechanism to explain the 



evolution of filament. Recent models of large scale corona (Mackay & van Ballegooijen 



2006a) consider the action of the large-scale surface motions such as differential rotation. 



meridional flow and supergranular surface diffusion. Our present observation indicates that 
large-scale surface flows are structured (not uniform) in the areas of divergence or stream 
flows, which should be taken into account in the numerical simulations. 

A better understanding of the mechanisms which lead to filament eruptions requires 
simultaneous multi-wavelength and multi-spatial resolution observations (both high-resolu- 
tion of the filament and low-resolution of the full Sun) over a wide range of latitudes. 
Indeed, previous works showed that different phenomena are observed at high resolution 



such as magnetic reconnection close to the starting location of the filament eruption (Rondi 



et al. , 2007). In this study, observing a larger area at lower resolution it was shown that 



at the same location where the filament first begins to erupt. We have found that a steep 
gradient in differential rotation, a north-south stream, and a shear in the zonal component 
were present. 
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8 Meridional magnetic flux transport 



8.1 Introduction 

The largest-scale velocity fields on the Sun consist of the differential rotation and meridional 
circulation. Both, the differential rotation and meridional circulation, are key ingredients 
of the solar dynamo. The differential rotation plays an important role in generating and 
strengthening toroidal magnetic field inside the Sun, while the meridional flow transports 
the magnetic flux towards the solar poles resulting in cyclic polar field reversals (for a 
recent review, 



see 



Brandenburg & Subramanian 2005 ) . 



The meridional flux transport seems to be an essential agent influencing the length, 
strength and other properties of the solar magnetic cycles. Generally, the slower the 
meridional flow, the longer the next magnetic cycle can be expected. Dynamo models 
showed that the turn-round time of the meridional cell is between 17 and 21 years, and 
that the global dynamo may have some kind of memory lasting longer than one cycle 



(Dikpati et al., 2006a). 



The speed of the meridional flow and its variation with the solar cycle measured by local 
helioseismology in the subsurface layers of the Sun are used as an input in the recent flux- 
transport models (Dikpati & Oilman 2006). In local helioseismology measurements (e.g. 



Zhao & Kosovichev 2004 Gonzalez Hernandez et al. 2006), the meridional flow is derived 



from a general subsurface flow field by averaging the north-south component of plasma 
velocity over longitude for a Carrington rotation period. The studies reveal that the mean 
meridional flow varies with the solar activity cycle. These variations may significantly 
affect solar-cycle predictions based on the solar dynamo models, which assume that the 



magnetic flux is transported with the mean meridional flow speed (Dikpati & Gilman 



2006). 



Our goal is to verify this assumption and to investigate the relationship between the 
subsurface meridional flows and the flux transport. In this study, we show that the mean 
meridional flows derived from the time-distance helioseismology subsurface flow maps are 
affected by strong local flows around active regions in the activity belt. However, these 
local flows have much less significant effect on the magnetic flux transport. 



8.2 Method of measurements 

The magnetic field data were obtained from Kitt Peak synoptic maps of the longitudinal 
magnetic field. The magnetic butterfly diagram is continuously constructed from synoptic 
magnetic maps measured at National Solar Observatory by averaging the magnetic flux in 
longitude at each latitude position for each solar rotation since 1976. 

At mid- latitudes of the magnetic butterfly diagram (Fig. 43 upper panel), between the 
active region zone and polar regions, we clearly see elongated structures corresponding to 



★ This chapter was submitted to Astrophysical Journal Letters as Svanda, M., Kosovichev, A. G., Zhao, J., 
Speed of meridional flows and magnetic flux transport on the Sun. 
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Figure 43: Upper panel - The magnetic butterfly diagram. Middle panel - The filtered magnetic 
butterfly diagram showing enhancements of the flux transport elongated structures. Bottom panel 
- The measured meridional flux transportation speed (in the south-north direction) for Carrington 
rotations 1900-2048. 



the poleward magnetic flux transport. The aim of our method is to measure the slopes of 
these structures and to derive the speed of the meridional magnetic flux transport. 

In addition to the large-scale structures, the original diagram contains small-scale rela- 
tively short-living local magnetic fleld structures, which appear as a 'noise' in the diagram. 
To improve the signal-to-noise ratio for the magnetic flux structures we applied a fre- 
quency band-pass fllter for the frequencies between 1.06 x 10~® s~^ (period of 1093 days) 
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Figure 44- Cartoon describing the method of measurement of the meridional flux transport speed 
from the magnetic butterfly diagram. 



and 3.17 x 10^^ s^^ (period of 36.5 days). The filtering procedure is performed sepa- 
rately for each individual latitudinal cut on the diagram. We tried also other methods of 
enhancement of structures and found that they all provide comparable results. 

The difference between the original butterfly diagram and the filtered one can be seen 
in Fig. 43 (upper and middle panel). The flux-transport elongated structures are more 
easily visible after the filtering, and therefore more suitable for analysis. 

The meridional flux transport is measured on the basis of cross-correlation of two 
latitude rows in the map. We assume that the fiux-transport structures in two different 
rows on the same hemisphere are similar in shape, but their positions are different due 
to the meridional transport. We cross-correlate pairs of rows separated by heliographic 
latitude Ab in a sliding window with the size of 55 Carrington rotations. The edges of the 
correlation window are apodized by a smooth function to avoid the boxcar effects. The 
extremal position is calculated as a maximum of the parabolic fit of the set of correlation 
coefficients of correlated windows in five discrete displacements. If the distribution of the 
correlation coefficients does not have a maximum, or if the normalized quality of the fit is 
too low (under 0.8), the meridional velocity in this pixel is not evaluated. We have chosen 
Ab = 5° as the best tradeoff between the spatial resolution and precision. 
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To make the procedure more robust, we average the calculated meridional velocity for 
five consecutive frames separated by 0.5° and centered at Ab from the studied row. If any 
of the speeds in averaged five rows is far out of the expected range (—60 to +60 ms~^), 
then it is not used in the averaging. From the fit, the accuracy of the measured fiow 
speed is evaluated and the maximum value of the set of five independent measurements at 
different rows is taken. 

The same procedure is done with the processed map rotated by 180° to avoid any 
possible preferences in the direction determination, and both results are averaged. The 
measured errors were taken as the maximum value of both independent measurements. 

For the ongoing analysis, only the speeds that were measured with the error lower than 
3 ms~^ were taken into account. This criterion and some failures of the slope measurement 
introduce gaps in the data, which we need to fill. For this purpose we need to determine 
the best continuous differentiable field that approximates the data. The determination of 
such a field can be done in various ways, but we wish to avoid possible artifacts. For filling 
the gaps we used the MultiResolution Analysis. It is based on wavelet analysis and we have 
chosen the Daubechies wavelet due to its compact support. This property is important 
since it minimizes edge effects. Moreover, using these wavelets also preserves the location 
of zero-crossing and maxima of the signal during the analysis. Daubechies are claimed to 
be very stable in the noisy environments. For details see Rieutord et al. (2007). 

For the comparison between the meridional flow obtained from time-distance helioseis- 
mology (Zhao & Kosovichev, 2004) and the magnetic flux transport from our method, we 
have calculated the averaged values of both quantities in bins of 10 heliographic degrees. 
For the period of 1996-2006, only eleven Carrington rotations have been evaluated (one 
per year) by time-distance helioseismology using the full-disc Dynamics data from the MDI 
instrument on SoHO spacecraft. These data are available only for 60 days per year. We 
compared the measurements of the magnetic flux transport and the meridional flows in 



those particular non-consecutive Carrington rotations. The plots are displayed in Fig. 45 



8.3 Results 

The correlation between the magnetic flux and flow speed is high (correlation coefficients 
are in range of 0.7-0.9). We have to keep in mind that while the time-distance meridional 
circulation proflles represent the behaviour of the plasma during particular Carrington 
rotations, the magnetic butterfly diagram tracking proflles represent the flux transport 
smoothed over 10 Carrington rotations. Therefore, the agreement cannot be perfect in 
principle. To make our results more accurate, the continuous helioseismic data are needed. 

The speeds of the meridional flux transport in the near-equatorial region are less re- 
liable, since the elongated structures in the magnetic butterfly diagram extend from the 
activity belts toward poles. In the equatorial region, signiflcant part of the measurements 
were excluded from the analysis due to their large measured error. The gaps were filled 
using MultiResolution Analysis from well-measured points. Although the results seem 
reasonable here, their lower reliability has to be kept in mind. 
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Figure 45: The longitudinally averaged meridional flow speed measured for a set of Carrington 
rotations by time-distance helioseismology. Solid line plots the time-distance mean meridional 
flow at 3-4-5 Mm depth, dotted at 6-9 Mm, and dashed at 9-12 Mm. The dots with error-bars 
represent the 10- degree-bin- averaged values of the flux transport speed derived from the magnetic 
butterfly diagram (Fig 43). 
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Figure 46: Left - A large-scale flow map at depth 3-4-5 Mm for Carrington rotation No. 1975 
(April 2001) and the corresponding MDI magnetogram in the grey-scale background. Large-scale 
flows towards the equator in the magnetic regions are visible around the large active region. Right 
- The longitudinally averaged meridional circulation profile for the same Carrington rotation. 
The southern hemisphere depicts almost no magnetic activity, so the meridional circulation pro- 
file obtained by averaging the time-distance flow map (solid line) almost fits the magnetic flux 
transport profile (points with error-bars) there, while on the northern hemisphere they differ. Af- 
ter masking the magnetic regions on the northern hemisphere, the recalculated profile (dotted line) 
tends to fit the butterfly tracking one also on the northern hemisphere. 



In the minimum of activity (such as CR 1911, 1923, or 2032), the profile of the merid- 
ional flux transport speed is very consistent with the mean longitudinally averaged profile 
of the meridional flow from helioseismology. The best agreement if found for the depth 
9-12 Mm. This suggests that the flux transport may be influenced by flows in the deeper 
layers. 

With increasing magnetic activity in the photosphere of the Sun, the gradient of the 
mean meridional circulation proflle derived from time-distance helioseismology becomes 



steeper. This is consistent with the results obtained by numerical simulations by Brun 



(2004). The simulations show that with increasing magnetic activity, the Maxwell stresses 
oppose the Reynolds stresses, causing an acceleration of the meridional circulation and 
deceleration of the rotation in low latitudes. Our measurements show that the variations 
of the slope of the mean meridional flux transport speed in latitute are lower with the 
progression of the solar cycle. 

When the large-active regions emerge in the activity belt, the flow towards equator 
is formed on the equatorial side of the magnetic regions (see example of the subsurface 
flow map in Fig. 46 left). This equatorward flow acts as a counter-cell of the meridional 
flow (present at the same longitudes as the corresponding magnetic region) and causes 
a descrease of mean meridional flow amplitude in the activity belt. This behaviour is 
noticed in all studied cases recorded during eleven non-consecutive solar rotations, for 
which the Dynamics data useful for helioseismic inversion exist. Therefore the formation 
of the apparent counter-cell seems to be a common property of all large active regions in 
depths 3-12 Mm. 

Flows in this counter-cell do not influence the magnetic flux transport, which can be 
demonstrated when the magnetic region is excluded from the synoptic map (Fig. [46] right). 
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Figure J^l: Cartoon of the mean meridional flow derived as the average of the north-south com- 
ponent of helioseismology flow maps. Flow in the periods without significant magnetic activity 
(solid line), in periods of increased magnetic activity (dashed line) and when the counter-cell in 
the active regions are present (dotted line). Measured meridional flux transport speed remains 
more-or-less similar to the solid curve. 



The calculated meridional circulation profile is then closer to the profile of the meridional 
flux transport speed derived from the magnetic butterfly diagram. 

8.4 Conclusions 

We have compared the measurements of meridional speed derived by two different tech- 
niques: by time-distance local helioseismology and by measuring the flux transport speed 
using the magnetic butterfly diagram. We have found that the results agree quite well in 
general, but they differ in the regions occupied by local magnetic fields. The detail flow 
maps from helioseismology show that this is due to the presence of meridional counter-cells 
at the equatorward side of magnetic regions, which influence the time-distance derived 
meridional flow profile, but does not influence the magnetic flux transport. 

We have studied eleven non-consecutive Carrington rotations covering one solar cycle. 
The effect of the local flows around active regions and especially on their equatorial side 



is noticed in all the studied cases (Fig. 47). Therefore, this behaviour seems to be a 



common property of the subsurface dynamics around active regions located in the activity 
belt. However, we have to keep in mind that both datasets are not directly comparable, 
since the time-distance flow maps represent the behaviour of flows during one Carrington 
rotation, while the butterfly diagram tracking procedure provide results averaged over few 
Carrington rotations. The more homogenous data for local helioseismology are needed to 
study this effect in more detail. 

The results show that the speed of the magnetic flux transport towards the solar poles 
may significantly deviate from the longitudinally averaged meridional fiow speed derived 
from local helioseismology measurements, which are affected by local circulation flows 
around active regions in the activity belt. Therefore using the mean meridional flow proflle 
from helioseismology may bias the results of the dynamo models. 
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9 Conclusions and perspectives 

During my Ph. D. studies I have developed the method suitable for the measurements of 
horizontal velocity fields in the solar photosphere, which is based on the local correlation 
tracking (LCT) method. The method consists of several separated steps covering the data 
processing procedure, including noise removal and coordinate transformations. The values 
of free parameters were adjusted using the synthetic data with known properties. They 
come from the SISOID numerical simulation, developed especially for this purpousc. From 
the comparison between the model velocity field and the calculated one the calibration 
curve was determined. At this point I have confirmed the fact that the LCT method 
underrepresents the magnitudes of the velocity field due to the spatial averaging during 
the calculation. 

1 have also compared the results of the proposed method with the results of time- 
distance local helioseismology, the method, which is considered the most powerfuU in the 
topic of velocity fields measurements. It is shown that both methods reasonably match, 
which is encouraging for both methods. 

The long-term behaviour of the surface flow field was studied. The generally accepted 
long-term properties were confirmed by the proposed method. The periodic analysis of the 
mean velocity components in the equatorial region did not show any conclusive results. 
We have detected few suspicious periods that might be related to the global Rossby waves 
pattern, but we cannot confirm this idea from the present material. 

When a slightly modified version of the method was applied to the data describing the 
temporal evolution of solar photosphere under the filament in the time before and after its 
eruption, it was found that the topology of the surface flows have changed signiflcantly. 
All the measured changes cause a stretch in the place, where the filament eruption started. 
The measurements support the theory of influence of the coronal magnetic held topology 
by the photospheric motions. 

The study of the meridional flux transport measured from the magnetic butterfly di- 
agram showed that the flux transport speed may signiflcantly differ from the meridional 
speed of the plasma measured by the time-distance local hehoseismology. The results sug- 
gest that in the dynamo models, the mean latitudinally averaged meridional flow speed 
coming from helioseismology cannot be used directly as the quantity describing the merid- 
ional flux transport. The longitudinal structure of the meridional flows seems to be very 
important. 

9.1 Perspectives 

The material obtained during the work related with this thesis provide a great perspective 
in the future investigations. In the futiu^c studies we would like to use it particularly for 
the investigation of the flows in active regions and of the coupling between the magnetic 
and velocity flelds in magnetic areas. Especially, we will work on possible conflrmation of 
the phenomenon of the dynamical disconnection of the sunspots from their magnetic roots, 
for which we believe that the current dataset is suitable. 
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By now, we have processed a three-day series of the great active region NOAA 9393, the 
largest active region of the current solar cycle, which are waiting for the detailed study, in 
which we would like to incorporate also the three-dimensional magnetic field reconstruction. 

The filament eruption, which was studied in the Section [7| is a perfect material for the 
non-potential modelling of the coronal magnetic field and its evolution in time. We expect 
that it will provide a lot of information about the particular coupling of the magnetic field 
in corona and the horizontal photospheric large-scale flows. 

The data coming from the magnetic butterfly diagram tracking analysis are suitable 
as an input for the flux transport dynamo models. We have started to cooperate with 
Dr. Mausumi Dikpati, an author of the well-behaving flux tranport dynamo model, which 
is able to reproduce the properties of last 12 solar cycles. We expect that the input of the 
real measured meridional flux tranport speed, which is essential, instead of that measured 
by local helioseimology will improve the results and the ability of the magnetic activity 
forecast. The longitudinal structure of the meridional flux transport speed and meridional 
flows need also a careful study. 
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10 Appendix - Flows on the stellar surfaces? 

Large-scale flows play and important role in the dynamics of the surface and sub-surface 
magnetic field. Their observation and measurements in the photosphere of the Sun bring 
information that can help theorists to model the processes involved in the solar dynamo. 
Assuming that the same physical principles like in our Sun take place also in the interior 
of other stars, it would be interesting also to obtain some information about the behaviour 
of the plasma in the photosphere of other stars. 

Recently, the progress in the Doppler imaging method ( Vogt et al. , 1987 or Rice et al. 



1989 ) allows the acquisition of time series of the images mapping the evolution and motions 



of surface features on the stellar surfaces (e. g. Kovari et al. , 2004). Although the resolution 
of surface maps is low and they contain lots of noise, they can be attributed to the presence 
of solar-like active regions caused by the presence of the magnetic field. The large-scale 
low resolution magnetic maps were used for the measurement of the surface flow on the 



Sun by the LCT method by Ambroz (2001a). Following this similarity, we may try to use 



the LCT method to track the surface features in the series of the Doppler imaging maps 
to derive the surface flow structure in the photospheres of other stars. 



10.1 LQ Hydrae 

As the pilot target the solar-like rapid rotator LQ Hydrae, for which a time-series of 
consecutive surface Doppler maps exists (Kovari et al. , 2004), was chosen. LQ Hya can 



be regarded as a very young proxy of our Sun. Astrophysical properties for LQ Hya are 
adopted from Kovari et al. (2004) and references therein. Relevant parameters for the 



present study are summarized in Table |3j 

Spot variability on LQ Hya has been studied since the discovery of its light variability. 



Strassmeier et al. (1997) and Strassmeier et al. (1999) detected rapid spot evolution on 



time-scales from weeks to months and estimated a possible spot cycle length of 6-7 years. 



With recent updates of the photometric data base to a total length of 21 years, Kovari 



et al. (2004) suggested an even shorter cycle length of 3.7 years, in agreement with the 



fundamental-mode oscillation period of 3.2 years predicted by Kitchatinov et al. (2000) 



from a distributed-dynamo model applied to LQ Hya. The first Doppler image of LQ Hya 



was presented by Strassmeier et al. ( 1993b ) and afterwards the star has become a favored 



target of different Doppler imaging campaigns until nowadays (the early history of Doppler 
imaging studies for LQ Hya was summarized by Strassmeier, [2002 ). The most recent 
Doppler imaging results were published in Kovari et al. (2004), where a time-series analysis 



was performed. From the cross-correlation of many consecutive Doppler images a weak 
solar-type surface differential rotation was recovered. 



For the time-series Doppler images, Kovari et al. (2004) used 52 spectra taken during 
a 57 nights long observing run at NSO in November— December 1996. The data covered 



★ This work was done together with Zsolt Kovari from Konkoly Observatory, Hungary, and Klaus G. Strass- 
meier, Astrophysikahsche Institut Potsdam, Germany, with contribution of Michi Weber, KataHn Olah, 
Kristian Vida, and Janus Bartus. 
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Parameter 


Value 


Classification 


IVZ V 


Distance (Hipparcos) 


io.oO ± U.OO pC 


Luminosity, Ij 


u./(U It u.uuy LiQ 




ft.u HI u.o 


^eff 


5070 ± 100 K 


V sini 


28.0 ± 1.0 kms-i 


Inclination, i 


65° ± 10° 


Rotation period, Prot 


1.60066 ±0.00013 days 


Radius, R 


0.97 ±0.07 i?0 


Age 


~ 100 Myr 


Mass 


~ 0.8 Mq 



Table 3: Astrophysical data for LQ Hydrae. Adopted from Kovdri et al. (2004) 0''nd Ferryman 



et al. (1991). 



35 consecutive rotations of LQ Hya and were used to recover a series of 28 overlapping 



Doppler images with the inversion code TempMap (Rice & Strassmeier, 2000 and Rice 



2002). The Fe 1-6430 A line is more reliable mapping line for the cool and relatively small 
f sini-star LQ Hya and resulted in maps with higher resolution compared to Cal-6439A. 
Therefore, the iron maps for the present pilot study were adopted. All maps show spot 
activity preferably at low latitudes, between —20° and ±50°, with a concentration in a 
band around the star centered at ±30°, and with only occasional evidence for a higher- 
latitude spot extension. Further information on the data, its reduction and preparations, 
and on the limits of the Doppler images of LQ Hya, can be found in the original paper by 



Kovari et al. (2004). 



The Doppler surface maps of LQ Hya were tracked with a fixed rotation period of 
1.60066 days (cf. Table [3]). Each of the 28 maps was transformed to the Sanson-Flamsteed 
coordinate system. The advantage of using sinusoidal projection is the conservation of 
areas. Only then maps are suitable for the application of LOT for the optimization of 
its best free parameters. Usually it is a trade-off between the sensitivity, i. e. the spatial 
resolution, and the signal-to-noise ratio - or equivalent - in the maps. In a parameter study, 
the trade-offs were determined that set the size of the correlation window FWHM and 
the time lag At between correlated images and found the best values with FWHM = 10° 
and At = 1 frame. These tests showed that the results are not very sensitive to the 
selection of FWHM but they are sensitive to the selection of At, which means that the 
difference structures are probably driven by rapid spot evolution not resolved in phase- 
averaged Doppler images. The mean time lag between consecutive images is 1.660 days, 
which equals to 1.037 rotations. Such consecutive images differ by only a single spectrum 
and differences between them are hardly representative. However, the trends of evolution 
are preserved in the entire time series. 
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Figure J^8: Top; Flows in the photosphere of the rapidly rotating K-dwarf LQ Hya (arrows). The 
length of the arrows is scaled by the actual radius of the stellar disc in solar units which is adopted 
to be O.97i?0. The background image is the segmented average temperature map used in the series 
of 28 Doppler maps to recover the flow field. We have truncated its temperature plotting range 
to emphasize the positions of its spots. There appear to be two distinct flow structures; a general 
poleward flow and a network of convergent flows around spots. Bottom; The root-mean- square 
deviations of individual flow vectors. 



10.1.1 Flow field 



For the entire series of 28 maps, 27 correlation pairs were calculated and averaged with 
equal weight. The resulting difference network is then interpreted as the surface flow field 
of LQ Hya and is shown as a vector field in Fig. |48] The average Doppler image of the 



series of 28 frames is used as the reference frame and is plotted as a background. The 
calculated root-mean-square error of individual vectors is plotted in the bottom panel of 
the same figure. The vectors in Fig. 48 indicate that the large-scale flow network is possibly 
related to the spatially resolved surface structures, which is expected from solar analogy. 
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Figure 49: Left; Mean azimuthal component (full line) of the flow field from Fig. ^8. The profile 
shows the signature of a solar-like differential rotation with higher angular velocity at the equator 
with respect to the poles. The best fit sin^ (p differential- rotation law is overplotted as a dashed 



line. Right; Mean meridional component of the flow field in Fig. 48 A northward flow prevailed 
on almost the entire stellar disc during the imaging period in 1996. As a comparison, the solar 
meridional flow is over-plotted as a dashed line. 

The amphtude of the velocity vectors of typically several hundreds of ms~^ also seems 
reasonable. Flows of the order of 100 ms~^ are observed in the photosphere of the Sun. 



Fig. |48] is indicative of the existence of convergent flows towards the starspot regions. 
Again, this is consistent with the situation in the solar photosphere, where convergent 
flows towards active regions are seen in the results obtained by local helioseismology (e. g. 



Zhao & Kosovichev, 2004). However, our usage of indirect Doppler maps of low temporal 
resolution incorporates that these convergent flows may not be uniquely interpretable as 
velocity fields and could also be due to the waxing or waning of spots. Because of the low 
spatial resolution of indirect stellar surface maps one cannot distinguish between a single 
large spot or a group of smaller spots. Consequently, we could have recorded a systematic 
migration pattern of small spots towards the center of an activity herd, e. g. in a giant 
convection cell. Nevertheless, it presents the first evidence and stimulates further study. 

10.1.2 Mean flows 



The average azimuthal velocity profile is displayed in Fig. |49] in the left panel. The curve 
can be fitted by a solar-like differential-rotation law of 

= fieq - AfisinV, (34) 

where VL{ip) is the angular velocity at a given latitude ip, VL^q is the equatorial rotation and 
Af2 is the rotational shear between the equator and the poles. The best fit is obtained 
with fieq = 3.928 ± 0.001 rad day"^ and Afi = 0.019 ± 0.002 rad day-\ in accordance 
with a lap time of ~ 330 days, i. e. the time the equator needs to lap the pole by one 
full rotation. As expected, the equatorial rotation rate is in perfect agreement with the 
phasing period of the Doppler images. The surface shear is consistent with that found by 
ACCORD (Average Cross-CORrelation of contiguous Doppler images) technique in the 
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study of |K6vari et aL] ( |2004[ ): Afi = 0.022 ± 0.008 rad day"^ After all, the two results 
from the two different methods confirm the principal reliability for both methods. 

The average meridional flow profile is displayed in Fig. 49 in the right panel. A poleward 
flow pattern seems to prevail for almost the entire stellar surface. The velocity peak near 
a latitude of +20° in Fig. 49 right is due to the poleward components of the converging 
flows around the many low-latitude spots. The underlying global meridional component 
may hkely just be seen above latitudes of, say, +40°. Fig. 49 left indicates a peak velocity 
of 0.01 rad/day (80 ms~^) at a latitude of above 60°. 

Rapid poleward migration of spots was first recognized by Vogt et al. (1999) on the 
sub-giant V711 Tau (HR 1099), and later verified through time-series Doppler imaging 
by Strassmeier & Bartus (2000). In both studies it was associated with the anti-solar 
differential rotation, i. e. the poles rotating faster than the equator. This was also seen on 
another giant, HD 31993 (Strassmeier et al. , 2003), and the issue was discussed recently 
by Weber (2006) suggesting that the stronger the anti-solar differential rotation, the faster 
the meridional flow towards the poles. On the Sun the meridional-fiow pattern is generally 
poleward and symmetric to the equator on both hemispheres. At that point, it has to 
be noted that Doppler reconstructions are more reliable near the visible polar hemisphere 
and become less reliable approaching the invisible part of the star. Therefore, detecting 
a solar-like north-south symmetry from Doppler images is beyond reason. On the other 
hand, variability of the differential rotation of active stars, including LQ Hya, is found on 
a timescale of few years by Petit (2006). 

Possibly connected, meridional flows on LQ Hya may also exhibit temporal fluctuations, 
just like on the Sun where the meridional flow is also observed to change its direction (e. g. 
Ulrich et al. , 1988). Therefore, the findings should be seen as a snapshot of an always 



changing pattern. More data in a further study will hopefully show more details on the 
general properties of the meridional flow of LQ Hya. 



10.1.3 Estimation of noise and sensitivity to free parameters 



Since the flow map in Fig. 48 is the average of 27 correlated pairs of frames, it is important 
to know the sensitivity to the selection of frames. The flow map from even- and odd- 
numbered frames was separately calculated, i. e. with time steps At of two frames. The 
flow maps from these two data sets are very similar, their median root-mean-square error 
for the whole map is only 10 ms~^. 

In the data set, two neighboring Doppler maps are not independent, they always differ 
by one spectrum (and overlap with usually seven spectra). As another test of reliability, 
the flow map using just the fully independent Doppler maps was calculated, i. e. running 
the calculations with At = 8 frames. In detail the flow map obtained this way is less 
reliable and more noisy. Although the surface features changed during the eight rotations 
covered, the median root-mean-square of difference of this map and the map calculated 
with At = 1 is still approximately only 50 ms~^. 
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10.1.4 Discussion 

This is a very preliminary study, which shows that the apphcation of LCT on the Doppler 
imaging data might work. It is clear from the root-mean-square errors displayed in Fig. |48] 
in the bottom part that the results are quite noisy. At this moment, the calculated flow 
field cannot be considered as reliable. There are a few points, at which the technique could 
be significantly improved. 

• The Doppler images used in the pilot study are reconstructed on the basis of only one 
spectral line. It would be more reliable to use during the inversions as many spectral 
lines as possible. The new version of TempMap code, which is in development, should 
enable this and therefore produce more reliable and less noisy temperature maps. 

• The used LCT code is the same as used for solar studies. Its use does not completely 
satisfy the assumption, under which it is considered reliable. For example the size 
of the correlation window is not larger than the tracers - starspots - but it is two 
or three times smaller. Although it has been tested on solar data that when using a 
sufficient averaging, also smaller FWHMs can be used, it is not clear enough, how 
this noise influences the calculated flow maps based on Doppler temperature maps. 

• The code does not take into account that temperature maps cover the whole stel- 
lar surface, not just the visible hemisphere. Therefore the whole code should be 
overwritten for this purpose. 

It is shown that an application of LCT technique on Doppler imaging maps provides mean 
flows comparable to those measured by different techniques. Therefore it can be concluded 
that the results are not completely arbitrary, however they may be arbitrary in details. This 
topic of research would be very important if it would be able to produce reliable surface 
flow maps that could be used as an observational input to the numerical simulations of 
global dynamo action. 



10.2 UZ Librae and a Geminorum 

UZ Lib is a KOIII RS CVn-type close binary variable, 3.4-times bigger than the Sun. 
The star is in detail described in the catalogue of chromospherically active binary stars 
by Strassmeier et al. (1993a). The period of rotation of the primary is the same as the 
orbital period of the system (Prot = -Porb = 4.768241 days) with vsini = 67 ± 1 kms~^, 
i = 50 °±10 °, and photospheric temperature Tes = 4800 K. In 1998, 11 consecutive Doppler 
images were obtained and now used for the recovering of the surface flows pattern. Each 6th 



-k The results from this section were presented at 5th Potsdam Thinkshop in June 2007 as parts of contri- 
butions K. Vida, Zs. Kovdri, M. Svanda, K. Oldh, J. Bartus and K. G. Strassmeier: Differential rotation 
and surface flow pattern on UZ Librae and Zs. Kovdri, J. Bartus, M. Svanda, K. Vida, K. G. Strassmeier, 
K. Oldh and E. Forgdcs-Dajka: Surface velocity network with anti-solar differential rotation on the active 
K-giant a Geminorum. 
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Figure 50: Surface flow pattern on UZ Librae, KOIII star. 
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Figure 51: Differential rotation and its parabolic fit and meridional circulation of UZ Lib 

and the solar meridional circulation (^rightj. 



frame was reconstructed from the independent set of spectral profiles in Fe 1-6393 A, Cal- 
6439 A, and Fe 1-6411 A lines; separately reconstructed images in each line were simply 
combined to get a better signal-to-noise ratio. The resulting flow map is displayed in 



Fig. [50] and a calculated differential rotation profile and meridional flow in Fig. [5T| For the 
meridional flow there are no convincing results, the error-bars are much larger than the 
mean value in each sampling point. The differential rotation is inconclusively anti-solar 
with parameters of fieq = 1-318 ± 0.002 rad day"^ and A^] = -0.005 ± 0.006 rad day-^ 
A stronger anti-solar rotation is found by ACCORD analysis of the same data set. 

For the calculation, just 5 maps were averaged, so that the signal-to-noise ratio is much 
lower than in the case of LQ Hya. Thus, the calculated maps are even less reliable, which 
depict itself in large errors of the differential rotation and meridional flow measurements. 

The star a Gem is a KlIII RS CVn-type long-period binary 12.3-times bigger than 
the Sun. As UZ Lib, a Gem is also locked (Prot = -Porb = 19.604 days) with t>sini = 
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Figure 52: Surface flow pattern on a Geminorum, KlIII star. 
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Figure 53: Differential rotation ('left ) and its parabolic fit and meridional circulation of a Gem 
( right;. 



27 ± 1 kms 



-1 



60° ± 15°, and photospheric temperature Teg = 4630 K. Together 6 



consecutive temperature maps were taken in 1997-1998 covering 3.6 rotations of the pri- 



mary component (Kovari et ah, 2001). Each second map was obtained using independent 



set of hne profiles in spectral lines Ca 1-6439 A and Fe 1-6430 A and combined together to 
decrease statistical uncertainties of individual Doppler reconstructions. 

Four independent pairs of frames were used for the LCT analysis and averaged to 



decrease the numerical noise in individual LCT maps. The flow map is displayed in Fig. 52 



The differential rotation (Fig. 53 left) shows anti-solar rotation with parameters of fieq = 

and Ai7 = —0.04 ±0.03 rad day~^. Almost ten times stronger anti 



13.83 ±0.01 rad day"^ and /\Vt = -0.04 ±0.03 rad day 
solar differential rotation is found using ACCORD analysis. The meridional flow (Fig. 53 
right) shows a poleward migration of the spot with an average speed of 220 ± 10 ms^^. 
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10.3 Conclusions 



The results in the field of stellar surface fiows investigation are in the state of individual 
pilot studies. There is a lot of issues, which has to be solved before the method will provide 
enough reliable results. With improving quality of the Doppler imaging procedures, for a 
lot of stars the series of evolving surface features may be expected in a short time period. 
The robotic telescopes such as STELLA (STELLar Activity; e. g. Strassmeier et al. , 2001 
Weber 2004 , or Bartus et al. 2006 ) will provide large amount of data suitable for stellar 



surface flow analysis. By that time it would be great to have a powerfuU tool for the stellar 
spots tracking and recovering the signal of differential rotation and meridional circulation. 
Recently, some theoretical dynamo models of main sequence stars including the differential 
rotation appeared (e. g. Lanza, 2006 or Cline et al. , 2003) and their results would be easily 



comparable to the measurements. 

Although the results obtained with the "standard" LCT code are not reliable, they are 
encouraging to develop the new code suitable for the tracking of features in the Doppler 
imaging temperature maps. Synthetic data would also provide important information 
about the reliability of the combination of Doppler imaging reconstruction and LCT ap- 
phcation. 

This topic of astrophysics will be continued in next years in cooperation with a group 
of Prof. Strassmeier in Potsdam. 
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